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1 .  INTRODUCTION 


It  is  now  generally  accepted  that  the  behaviour  of  a  reactive  material 
needs  to  be  considered  in  terms  of  its  explosive,  its  thermal  and  its  mech¬ 
anical  properties.  Most  explosives  deform  before  initiation  and  the  details 
of  this  deformatiou  determine  how  energy  can  be  localised  to  give  "hot  spots". 
Once  heat  is  produced  locally  the  "hot  spot"  temperature  is  controlled  by 
the  heat  being  produced  and  the  flow  of  heat  away  from  the  "hot  spot".  To 
understand  this  process  it  is  important  to  have  the  reaction  parameters  of 
the  explosive  and  thermal  properties  such  as  specific  heat  and  conductivity. 

For  many  explosives,  particularly  new  compositions,  these  are  not  always 
known  accurately. 

The  research  on  this  contract  has  considered  the  mechanical  and  thermal 
properties  in  some  detail.  Studies  of  initiation  and  propagation  have  been 
greatly  aided  by  high-speed  photography.  Research  on  reaction  kinetics  and 
thermal  properties  has  used  mass  spectroscopy,  differential  scanning  calori¬ 
metry  (DSC),  thermogravimetric  analysis  (TG)  and  scanning  electron  microscopy. 
Mechanical  properties  experiments  include  measurements  of  hardness,  yield 
strength,  coefficient  of  friction  and  fracture  surface  energy. 

The  four  areas  of  work  described  in  this  report  have  either  been 
published  (sections  2,3,4)  or  are  ready  for  submission  (section  I). 

The  texts  of  these  papers  are  given  in  full. 

2.  THE  IGNITION  OF  A  THIN  LAVER  OF  EXPLOSIVE  BY  IMPACT;  THE  EFFECT  OF  ADDED 
POLYMER  PARTICLES.  G.M.  Swallowe  and  J.E.  Field. 

2 . I  Introduction 

A  study  of  the  deformation  behaviour  of  a  thin  layer  of  material  when 
impacted  is  relevant  to  a  range  of  problems  including  the  sensitiveness  of 
explosives.  A  standard  test  procedure  for  assessing  the  hazard  involved  in 
handling  explosive  materials  is  to  impact  a  sample  (typically  20-50  mg)  with 
a  falling  weight.  The  sensitiveness  of  a  sample  is  usually  expressed  in 
terms  of  a  50%  drop  height  to  cause  ignition  (i.e.  the  height  from  which  the 
weight  would  ignite  50%  ol  a  series  of  samples).  There  are  standard  methods 
for  obtaining  the  50%  height  (sec,  for  example,  Dixon  and  Massey,  1957).  However 
although  the  test  has  been  used  for  many  decades  it  is  only  recently  that 
a  systematic  attempt  has  been  made  to  photograph  the  sample  during  its 
deformation  (Heavens,  1973,  Heavens  and  Field  1974).  A  combination  of 
high-speed  photography  combined  with  pressure  measuring  techniques  helped 
establish  the  physical  processes  occurring  during  ignition  and  propagation. 

In  general  it  is  thought  that  the  initiation  of  an  explosive  by  mech¬ 
anical  shock  is  thermal  in  origin,  although  some  workers  believe  that  a 
tribochemieal  or  a  moLecuLar  fracture  me eh  inism  may  be  responsible  in  certain 
circumstances  (Taylor  and  Weaie,  1932,  Ubbolohde ,  1948,  Fox,  1970).  On  the 
basis  of  localized  thermal  energy  or  "lot  spots"  as  the  source  of  the  ex¬ 
plosion,  four  possible  mechanisms  have  been  envisaged  for  ignition  by  impact. 

(i)  Adiabatic  compression  of  trapped  gas  spaces  (Bowden,  Mulcahy,  Vines 
and  Yoffe,  1947,  Chaudhri  and  Field,  1974 >. 

(ii)  Viscous  heating  of  material  rapidly  extruded  between  the  impacting 
surfaces  (Eirich  and  Tabor  1947,  Bolkhoviwinov  and  Pokhil  1958)  or  by  capi¬ 
llary  flow  between  grains  (Rideal  and  Robertson,  1948). 
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(iii)  Friction  between  the  impacting  surfaces  and/or  grit  particles 
and/or  grains  of  the  material  (Bowden  and  Gurton,  1949). 

(iv)  Localised  adiabatic  deformation  of  the  thin  layer  upon  mechanical 
failure  (Af anas' ev  and  Bobolev, 1971 ,  Winter  and  Field,  1975). 

It  was  shown  quite  convincingly  (Bowden  et  al,  1947)  that  the  presence 
of  gas  bubbles  affected  the  impact  sensitiveness  of  nitroglycerine, and 
Bowden  maintained  repeatedly  (Bowden  and  Yoffe,  1949,  1952,  1958,  Bowden, 
1950,  1963)  that  mechanism  (i)  was  the  primary  cause  of  ignition,  not  only 
in  liquids  but  also  in  solid  explosives.  This  view  has  with  equal  persis¬ 
tence  been  rejected  by  the  Russian  school  of  thought  (Andreev,  Maurina  and 
Rusakova  1955,  Bolkhovitinov,  1959,  Afanas'ev  and  Bobolev, 1971 ) .  In  impact 
experiments  on  solid  explosive  materials  it  was  discovered  (Kholevo,  1946) 
that  the  sensitiveness  was  considerably  reduced  if  the  sample  was  prevented 
from  flowing.  Afanas'ev  and  Bobolev (1961)  produced  evidence  from  strain- 
gauge  experiments  which  suggested  that  during  impact  the  sample  suddenly 
undergoes  a  type  of  fracture.  These  authors  proposed  a  mechanism  of  hot¬ 
spot  formation  by  the  release  of  energy  along  slip  surfaces  formed  at  the 
onset  of  mechanical  failure  (Afanas'ev  and  Bobolev* 1971 ,  Afanas'ev,  Bobolev, 
Kazarova  and  Karabanov  1972). 

The  high-speed  camera  work  of  Heavens  and  Field  (1974)  showed  that 
impacted  samples  in  the  drop-weight  test  may  undergo  bulk  plastic  flow, 
show  evidence  of  partial  fusion  and  even  (with  PETN)  melt  completely.  The 
flow  speed  during  these  processes  is  considerable  and  may  reach  a  few  100  m  s 
Ignition  occurred  at  a  small  number  of  local  hot-spots.  Strain  gauge  measure 
ments  showed  that  pressures  of  typically  0.5  -  1  GPa  (5-10  kbar)  were  ach¬ 
ieved.  If  a  sample  failed  by  plastic  flow  this  was  accompanied  by  a  sharp 
pressure  drop.  It  was  after  the  pressure  drop  that  the  high  velocity  flow 
phenomena  took  place  and  ignition  occurred.  Hot-spots  were  thought  to  be 
caused  either  by  local  obstructions  in  the  flow  or  by  the  rapid  closure  of 
gas  pockets  trapped  in  the  material. 


Recent  work  has  been  concerned  with  studying  the  effect  of  added  grit 
particles  in  greater  detail.  The  case  of  the  hard,  high  melting  point 
particle  is  reasonably  well  understood  (Bowden  and  Gurton  1949,  Bowden  and 
Yoffe  1952,  Ubbelholde  1948).  The  basic  idea  is  that  when  two  solids  rub  or 
impact  together  the  hot-spot  temperature  at  the  interface  is  determined  by 
the  solid  with  the  lowest  melting  point.  Properties  such  as  thermal  con¬ 
ductivity  and  hardness  are  important  but  only  as  second-order  effects.  In 
fact  by  chosing  grits  of  different  melting  point  and  measuring  the  impact 
sensitiveness  Bowden  and  Gurton  were  able  to  measure  hot-spot  ignition 
temperatures  for  a  range  of  explosives.  The  hot-spot  temperature  found 
by  these  workers  for  PETN  in  a  drop-weight  impact  situation  was  ca.  700  K 
for  hot-spots  of  micron  dimensions. 

It  is  now  known  that  certain  polymers  such  as  polycarbonate  CPC)  and 
polysulphone  (PS)  can  nlso  sensitise  explosives  in  the  drop-weight  test 
though  others  such  as  polypropylene  (I’P)  do  not  (Bean  1973).  At  first 
sight  this  result  appears  surprising  since  polymers  arc  relatively  soft 
and  have  softening  temperatures  usually  well  below  the  700  K  mentioned 
above.  The  primary  object  of  our  new  work  was  to  try  to  understand  the 
mechanisms  by  which  some  polymers  give  increased  sensitiveness  and  why  others 
do  not.  At  the  start  it  appeared  possible  that  the  explanation  could  be  a 
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chemical  one  based  on  the  production,  of  free  radicals  which  aided  the 
reaction  kinetics  of  the  explosive  or  a  mechanical  one  in  which  hot-spots 
were  produced  by  the  deformation  of  the  polymer.  The  first  experiments 
were  concerned  with  possible  chemical  effects  and  used  techniques  such  as 
differential  scanning  calorimetry  (DSC)  and  thermogravimetric  analysis  (TG) . 
The  mechanical  effects  were  studied  using  the  photographic  and  pressure 
measuring  techniques  described  earlier. 

2.2  Studies  of  Chemical  Effects 

Experimental 

These  experiments  involved  measuring  the  reaction  parameters  for  decom¬ 
position  of  pure  explosive  and  explosive/polymer  mixtures  and  assessing 
any  differences.  The  apparatus  used  was  a  Stanton-Redcroft  TG-750  thermo- 
balance  which  allows  the  weight  of  the  sample  to  be  monitored  whilst  it  is 
heated  at  a  constant  rate,  thus  giving  a  weight-loss  versus  temperature 
curve.  The  samples  used  consisted  of'M  mg  lots  of  a  mixture  containing 
25Z  (by  weight)  of  polymer  to  75X  explosive.  This  percentage  of  polymer 
is  much  greater  than  likely  to  be  found  in  practice,  but  was  chosen  to 
amplify  any  interaction.  The  explosive  used  in  all  these  experiments  was 
pentaerythritol  tetranitrate  (PETN).  The  analysis  follows  the  method 
of  Hauser  and  Field  (1978  and  section  3  of  this  report)  and  so  ia  only 
briefly  outlined  here. 

Decomposition  is  assumed  to  follow  an  equation  of  the  form 


—  -  -  A  exp(-E/KT)wn  (I) 

with  w  the  fractional  residual  weight  at  time  t,  E  the  activation  energy, 

A  the  frequency  factor,  n  the  reaction  order  and  T  the  temperature.  Since 
the  heating  rate  H  is  constant,  one  may  write  dT/dt  ■  H.  Substituting 
in  equation  1  and  taking  logs  yields: 

lnlf  *  lr4  "  It  +  n  ln  w  (2) 

Thus  a  plot  of  (ln  dv/dT  -  n  ln  w)  against  l/T  gives  a  straight  line  if 
n  has  been  chosen  correctly.  The  value  of  E  is  determined  from  the  slope 
of  the  line  and  A  from  its  intercept.  The  analysis  was  carried  out  on  the 
Cambridge  IBM  370/165  computer.  The  programs  used  fit  a  cubic  spline  to 
the  original  data  and  produce  differentials  which  are  then  used  to  give  plots 
of  (ln4s  -  n  in  w)  against  l/T  for  values  of  n  from  0  to  2  in  steps  of  0.2. 
The  straightest  section  of  these  curves  for  each  stage  of  the  reaction 
there  may  be  more  than  one)  was  then  selected  by  eye  and  the  appropriate 
values  of  n,  E  and  A  obtained.  The  method  is  illustrated  by  the  activation 
energy  plot  for  a  sample  of  pentaerythritol  tetranitrate  (PETN)  in  Figure  I . 
The  technique  was  assessed  in  two  ways.  Firstly  by  using  the  obtained 
reaction  parameters  and  simulating  a  decomposition  curve  which  could  be 
competed  with  the  original.  Secondly,  the  dehydration  of  calcium  oxalate 
monohydrate  was  observed. 

CaC20AH20  ♦  CaC204  ♦  H2<> 
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This  reaction  is  well  understood  and  the  reaction  parameters  are  well  docu¬ 
mented  in  the  literature  (Gurrieri,  Siracusa,  Cali,  1974,  Freeman,  Carroll, 
1958).  The  results  obtained  with  the  technique  outlined  above  were  in  good 
agreement  with  the  literature  values  and  are  set  out  in  Table  I . 

Table  1 ;  Dehydration  of  Calcium  Oxalate  Monohydrate 

E/kJ  mol  1  In  A/ln  e  \ 

This  work  89.5  -  3  20,2 

Literature  88.6  21.1 


Results 

(a)  Thermal  decomposition  of  PETN/polymer  samples 

The  polymers  reported  to  sensitise  PETN  (Bean,  1973)  were  polycarbon¬ 
ate  (PC),  polysulphone  amber  (PS)  and  to  a  lesser  extent  polypropylene  (PP). 
Results  obtained  in  this  work  confirmed  the  dramatic  effect  produced  by  PS 
and  PC  but  showed  very  little  increased  sensitivity  in  the  presence  of  PP  . 
Samples  of  these  polymers  supplied  by  AWRE  were  ground  in  a  freezer  mill  to 
sizes  of  <76  urn  and  then  mixed  with  PETN  to  produce  experimental  samples. 
The  heating  rate  used  throughout  was  10  deg  min  and  the  purger  gas 
nitrogen  or  argon.  The  results  are  sumnarised  in  Table  2,  the  figures 
listed  being  averages  of  between  three  and  seven  different  experiments. 


Table  2:  Effect  of  Polymeric  Additives  on  the  Thermal 
Decomposition  of  PETN 

Sample  E  /kJ  mol  1  E./kJ  mol""1  In  A  /In  s  1  In  A. /In  s  1  n  T/K 
s  d  s  d _ 


PETN 

143 

♦ 

5 

193  -  7 

31 

+ 

4 

44 

-  6 

0.7 

478 

PETN+PC 

130 

♦ 

10 

214  -30 

30 

+ 

6 

53 

-io 

1.0 

476 

PETN+PP 

122 

♦ 

15 

172  -15 

28 

+ 

4 

41 

-  6 

0.6 

477 

PETN+PS 

130 

♦ 

10 

231  -20 

30 

+ 

4 

49 

i  ♦ 

© 

1.0 

474 

In  each  case,  a  two  stage  reaction  was  observed  (Figure  1),  the  first 
stage  being  identified  with  the  sublimation  of  PETN  (E  ,A  )  and  the  second 
its  decomposition  (E,  ,A.).  The  values  obtained  for  the  aStivation  energies 
of  sublimation  043  *a5  KJ  mol-,)and  decomposition  (193  *  7  kJ  mol"!)  are  in 
good  agreement  with  the  literature  values.  Cundall,  Palmer  and  Wood, (1977) 
quote  150.4  kJ  mol”*  for  the  sublimation  of  PETN  while  Roberston  (1948), 
Rogers  and  Morris  (1966)  Maycock  and  Vemeker(l970)  and  Ng,  Field  and  Hauser 
(1976)  give  values  for  the  decomposition  of  197,198,188  and  192  kJ  aol“* 
respectively.  The  value  for  In  A^  of  43.5  is  also  in  good  agreement  with 
the  45.6  reported  by  Robertson  (1948). 

A  larger  spread  is  observed  in  the  values  of  in  the  PETN/additive 
samples  than  the  pure  PETN.  However,  although  some  activation  energy 
differences  are  observed,  there  are  no  significant  changes  or  consistent 
trends.  The  values  of  frequency  factor  are  much  less  accurate  than  those 
of  activation  energy  since  they  have  been  obtained  by  measuring  the  intercept 


*  The  sensitising  effect  reported  by  Bean  for  PP  in  PETN  is  surprising 
and  we  query  the  effect.  It  is  contrary  to  our  results  and  in  fact  Bean's 
results  for  PP  in  other  explosives  show  little  or  no  positive  effect.. 
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rather  than  the  slope.  Again, within  the  errors  of  measurement, they  show  no 
significant  changes  from  the  values  for  pure  PETN.  If  the  additive  had 
significantly  affected  the  reaction  towards  greatly  increased  reaction  rate 
then  A  should  have  increased  and  E  decreased.  Table  2  shows  that  only  small 
changes  in  A  and  E  took  place ,  and  further  that  any  changes  of  A  and  E 
tended  to  cancel  (i.e.  as  A  increases,  E  also  increases  and  vice  versa). 

The  only  consistently  observable  feature  in  the  results  is  a  change 
in  the  value  of  n  from  0.7  to  1.0  (except  in  the  case  of  PP).  A  value  of 
n  ca.  0.7  implies  a  surface  reaction  (for  which  n  =  2/3)  and  n  *  1 .0 
is  the  value  one  would  expect  for  a  bulk  reaction.  The  explanation  for 
this  is  that  the  high  proportion  (25%)  of  filler  material  distributed 
throughout  the  sample  breaks  up  the  PETN  giving  it  a  greater  surface  area, 
thus  favouring  a  bulk  rather  than  a  purely  surface  reaction.  The  result 
for  PP  in  explained  by  the  fact  that  it  softens  and  melts  at  a  lower  temp¬ 
erature  than  the  other  polymers.  PETN/PP  samples  were  observed  to  give 
droplets  of  polymer  which  then  coalesced  to  form  one.  drop.  This  effectively 
removed  the  dispersion  of  particles  which  gives  n  approaching  unity  leaving 
a  surface  reaction  with  n  ca.  2/3. 

(b)  Thermal  decomposition  of  PETN/Benzoyl  peroxide  samples 

The  results  of  the  preceeding  section  show  that  the  additives  tested 
have  very  little  effect  on  the  thermal  decomposition  of  PETN.  It  may, 
therefore  be  concluded  that  their  sensitising  effect  is  not  due  to  polymer 
decomposition  products  interacting  with  PETN.  However,  the  possibility 
still  exists  that,  when  polymers  fracture  or  flow  during  impact,  large  con¬ 
centrations  of  free  radicals  may  be  produced,  and  these  could  react  with 
PETN  and  cause  ignition.  Walker  and  Green  (1976)  have  proposed  that  the 
mechanical  production  of  free  radicals  is  important  in  initiation  and  have 
done  some  experiments  in  which  they  find  that  tetraroethylamroonium  tribo- 
hydride  (a  free  radical  donor)  causes  increased  reaction  with  ammonium 
nitrate. 


In  order  to  test  this  hypothesis,  mixtures  of  benzoyl  peroxide  and 
PETN  were  used  in  a  series  of  TG  experiments.  Benzoyl  peroxide  is  an 
organic  oxide  which  thermally  breaks  down  in  the  temperature  range  310 
to  350  K  to  form  two  identical  free  radicals 
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It  is  widely  used  in  the  polymer  industry  as  a  ’starter’  material  for 
free  radical  polymerisation  as  it  provides  a  good  source  of  easily  formed 
free  radicals. 


A  number  of  TG  runs  were  carried  out  on  both  benzoyl  peroxide  by 
itself  and  two  PETN/benzoyl  peroxide  mixtures.  The  mixtures  used  had  com¬ 
positions  of  75%  to  25%  and  93%  to  7%  by  weight  of  PETN/benzoyl  peroxide. 

Due  to  the  decomposition  of  the  benzoyl  peroxide,  the  only  activation 
energy  that  can  be  estimated  is  the  decomposition  of  PETN,  the  sublimation 
stage  being  buried  in  the  benzoyl  peroxide  reaction.  Simulated  computer 
decompositions  of  these  mixtures  were  made  by  combining  75%  of  the  numerical 
value  of  pure  PETN  run  with  25%  of  the  numericr 1  value  of  a  pure  benzoyl 
peroxide  run  (and  similarly  for  the  93/7  mix}.  The  figures  thus  obtained 
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were  fronted  ns  if  they  had  been  real  (lain  and  the  analysis  carried  out  in 
the  usual  manner.  Comparison  of  the  K,A  and  n  values  obtained  in  real 
and  simulated  experiments.  Table  3,  reveals  that  what  appears  at  first  to 
be  a  very  strong  chemical  interaction  is  due  to  the  combined  effects  of 
the  two  materials  acting  independently. 
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(c)  Effect  of  ultra-violet  light  on  the  thermal  decomposition  of  PETN /grit 
samples 

As  a  further  test  of  the.  free-radical  initiation  hypothesis,  PETN/ 
additive  samples  were  irradiated  with  ultra-violet  light  as  they  v?ere  being 
heated  in  the  TC.  ll-V  is  frequently  effective  in  degrading  polymers  by 
breaking  chains  with  the  formation  of  free  radicals.  Both  PP  and  PS  arc 
susceptible  to  photodegrad.it  ion  by  If— V  whereas  PC  is  degraded  relatively 
much  less.  It  would  therefore  be  expected  that,  if  free  radicals  are 
important  in  sensitising  PETN ,  both  PP  and  PS  in  combination  with  U-V  should 
show  an  appreciable  effect  and  PC  little  or  no  effect.  The  results  obtained 
.in  these  experiments  are  set  out  in  Tabic  4. 
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Except  for  the  rather  anomalous  behaviour  •  >!  PS  in  the  decomposition  stage 
(see  later),  the  results  indicate  that  •>n)  ’  small  effects  ate  produced. 


The  changes  in  the  sublimation  parameters  E  and  A  tend  to  cancel  in  all 
cases  and  for  the  decomposition  parameters  and  changes  again  effect¬ 
ively  compensate,  particularly  for  polycarbonate.  PP,  PS  and  pure  PETN 

all  show  decreases  in  (E^  in  the  case  of  PS)  combined  with  increases 
in  (Ad2  with  PS).  Thus  the  reaction  is  helped,  albeit  a  little,  by 
the  U-V  radiation.  This  is  illustrated  by  the  slight  but  consistent, 
reduction  in  the  reaction  peak  temperatures  Tp. 

PS  shows  the  most  significant  changes  in  that  an  extra  reaction  stage 
appears.  This  result  was  found  to  be  very  repeatable,  though  it  is  not 
fully  understood.  Although  Ed |  (for  the  new  stage)  is  very  high,  the  value 
of  Ad |  more  than  compensates  for  this,  and  Tp  is  reduced  by  12  K.  The 
initial  and  final  stages  of  the  reaction  are  virtually  unchanged  by  the 
presence  of  this  new  stage,  which  only  has  a  small  (ca.  9  K)  range  of  activity. 

A  further  series  of  experiments  on  the  PETN/PS  system  was  carried  out 
in  which  the  U-V  light  was  switched  on  or  off  during  the  course  of  a  run 
and  the  effect  on  the  reaction  parameters  noted.  These  experiments  show 
that  the  stage  of  reaction  at  which  the  U-V  is  present  is  quite  important. 

If  the  U-V  is  switched  on  at  any  stage  before  a  third  of  the  material  present 
has  decomposed,  then  a  three  stage  reaction  is  observed,  otherwise  no  effect 
is  observed.  When  U-V  has  been  present  for  the  initial  stages  of  the  reaction 
(even  if  only  one  eighth  of  the  explosive  has  decomposed)  extinguishing  the 
light  has  no  effect  and  the  reaction  proceeds  in  three  stages,  as  if  the 
U-V  were  present  throughout. 

These  results  indicate  that  the  U-V  is  effective  in  breaking  down  the 
PS  to  produce  free  radicals  which  then  influence  the  decomposition  of  PETN. 

If,  however,  PETN  decomposition  has  proceeded  to  a  significant  extent  before 
the  U-V  is  present,  the  irradiation  has  little  effect.  This  is  probably 
because  the  free  radical  concentration  produced  by  the  PETN  alone  is  already 
much  greater  than  that  which  the  lamp  caiproduce.  The  observation  that 
extinguishing  the  lamp  early  in  the  reaction  does  not  stop  the  three  stage 
process  indicates  that  the  major  influence  is  the  concentration  of  free 
radicals  in  the  initial  stages  of  decomposition.  The  influence  of  U-V  may 
thus  be  summarised  as  producing  such  a  concentration  of  radicals  in  the 
early  stages  of  PETN  decomposition  that  the  decomposition  process  is  speeded 
up  and  the  reaction  parameters  (and  presumably  the  details  of  the  mechanism) 
are  changed. 

That  PC  is  not  influenced  in  this  way  is  not  surprising  since,  in 
comparison  with  other  polymers,  polycarbonates  are  not  degraded  to  any 
large  extent  by  U-V  (Ranby  and  Rabak,  1975).  It  is,  however,  surprising 
that  if  such  an  effect  should  be  produced  with  PS,PP  remains  ineffecitve. 

Both  PS  and  PP  are  degraded  to  a  large  extent  by  U-V  and  both  produce 
large  concentrations  of  radicals  when  degraded.  It  suggests  that  the 
sulphone  groups  have  a  specific  effect  on  at  least  one  of  the  stages  of  the  PETN 
decomposition.  Support  for  this  comes  from  the  work  of  Reich  (1973)  who  tested 
the  compatibil ity  of  a  range  of  polymers  with  various  explosives.  PS  is  one 
of  the  polymers  he  lists  as  incompatible  with  PETN  since  he  detected  an 
increase  in  the  value  of  the  reaction  parameter  n  in  the  present  of  PS. 


From  the  series  of  experiments  described  in  this  section,  it  may  be 
concluded  that  there  is  no  major  chemical  reaction  between  PETN  and  any 
of  the  polymer  fillers.  The  only  exception  to  this  is  PS  and  then  only  in 

the  extreme  conditions  of  having  a  high  intensity  of  U-V  shone  on  the  sample. 
However,  the  effect,  even  with  PS,  is  still  only  marginal  and  the  conclusion 
has  to  be  that  chemical  effects,  although  present,  are  only  of  secondary 
importance  compared  with  the  mechanical  ones  to  be  discussed  in  the  following 
section . 


2 .  1  Photographic  Observations  of  the  Impact  Behaviour  of  Polymers 
and  Explosives 

(a)  High-speed  photography  of  impacts  on  polymer  samples 

The  impact  behaviour  of  both  polymers  and  explosives  was  photographed 
using  a  C4  rotating  mirror  framing  camera  at  framing  intervals  of 

ca.  6.5  ps.  Impacts  were  observed  using  a  transparent  toughened  glass 
anvil  system  similar  to  that  of  Heavens  and  Field  (1974)  and  is  shown 
schematically  in  Figure  2.  The  weight  of  5.5  kg  was  dropped  by  an  elec¬ 
tromagnet  from  a  height  of  1  m,  the  fall  being  guided  by  three  cylindrical 
rods . 


The  behaviour  of  small  discs  of  polymer  ca.  2  mm  diameter  and  0.8  mm 
height  under  impact  is  illustrated  in  Figure  3.  Figure  3  shows  selected 
frames  from  an  impact  on  a  disc  of  PP.  A  trace  of  radial  expansion  versus 
time  is  given  in  Figure  5.  With  this  material,  there  are  no  rapid  changes 
in  the  expansion  rate  and  bulk  deformation  occurs  throughout.  Figure  4  is 
a  higher  magnification  sequence  for  an  initially  2  mm  diameter  disc  of  PS. 
Its  radial  expansion  is  also  included  in  Figure  5.  This  material  deforms 
initially  in  a  bulk  manner  but  eventually  fails  catastrophically  (see 
frame  4(c)  and  Figure  5).  PC  gives  similar  catastrophic  failure  to  PS. 
Microscopic  examination  of  the  samples  after  impact  confirms  that  catas¬ 
trophic  failure  of  PC  and  PS  is  associated  with  rapid  cracking  and  shearing 
in  localised  bands.  Figure  6  is  a  micrograph  of  a  deformed  PS  disc.  There 
is  a  network  of  fracture  throughout  the  sample  and  sets  of  fine  parallel 
shear  bands  (one  set  is  labelled  s).  An  enlarged  view  of  region  s,  but 
rotated  so  that  the  bands  are  horizontal  in  the  figure  is  given  in  6b.  It 
is  argued  later  that  it  is  these  regions  of  localised  deformation  (crack¬ 
ing  and  shear  banding)  which  give  rise  to  "hot-spots”. 

These  observations  lead  to  the  conclusion  that  the  polymers  used  may 
be  classified  in  two  broad  groups:  (i)  Polymers  such  as  PC  and  PS  which  fail 
by  the  production  of  many  fast  moving  cracks  and  local  shear  bands.  This 
is  accompanied  by  rapid  radial  expansion.  (ii)  Polymers  such  as  PP  which 
deform  plastically  at  quite  a  high  rate  but  do  so  by  bulk  deformation  and 
without  cracking  and/or  shear  banding. 

(b)  High-speed  photography  of  impacts  on  explosives 

Heavens  and  Field  (1974)  photographed  the  impact  of  a  wide  range 
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of  granular  explosives.  Powdered  PETN  was  observed  to  be  compacted  into 
a  pellet  and  then  to  undergo  severe  plastic  deformation.  The  material 
became  translucent  and  finally  completely  transparent  with  very  rapdid  f low 
(300  m  s~l)  which  Heavens  and  Field  attributed  to  melting.  Ignition  occurred 
in  this  final  stage.  RDX  partially  fused  when  impacted  between  glass  anvils 
and  there  was  indirect  evidence  of  complete  sample  melting  of  RDX  between 
steel  anvils.  Pressure-time  curves  obtained  by  Heavens  and  Field  showed 
that  a  sharp  pressure  drop  occurred  when  the  sample  failed  plastically. 
Ignition,  in  samples  without  added  grit  particles,  always  occurred  after 
the  pressure  drop  and  when  the  explosive  was  flowing  rapidly.  In  the 
present  work  a  large  number  of  experiments  were  carried  out  on  pellets  of 
PETN  and  the  behaviour  was  very  similar  to  that  described  by  Heavens 
and  Field  for  granular  material.  Figure  7  illustrates  a  typical 
result  in  which  a  25  mm  pellet  of  PETN,  pressed  to  5k  bar  and  of  diameter 
5  mm  is  impacted  by  a  5.5  kg  mass  from  a  height  of  1  m.  The  same  type  of 
behaviour  as  described  for  powders,  can  clearly  be  observed.  In  frames 
a  to  f  the  material  is  deformed  plastically  and  becomes  translucent.  Trans¬ 
parency  and  rapid  flow  (230  ms’)  occur  in  frame  g  with  ignition  at  sites  A 
and  B  almost  simultaneously. 

The  results  obtained  with  a  large  number  of  such  sequences  show  that 
ignition  only  occurs  if  rapid  flow,  as  observed  in  frame  g,  takes  place. 
Although  rapid  flow  is  a  necessary  prerequisite  for  ignition  the  converse 
is  not  necessarily  true  and  many  failures  occur  without  a  resulting  ignition 
This  suggests  that  there  has  to  be  a  particular  region  in  the  flowing  material 
where  a  "hot-spot"  develops.  These  regions  could  be  associated  with  gas 
bubbles,  foreign  particles  or  enhanced  flow.  If  discontinuities  of  sufficient 
size  do  not  exist,  ignition  does  not  take  place.  Note  that  fast  reaction 
always  starts  at  localised  regions  (see  Figure  7  and  the  sequences  in 
Heavens  and  Field). 

( c )  High-speed  photography  of  impacts  on  explosive/polymer  samples 

A  number  of  high-speed  photographic  sequences  at  framing  rates  of  up 
to  2  x  10^  frames  per  second  have  been  taken  of  impacts  on  explosive/grit 
samples.  These  sequences  were  taken  in  transmitted  light  using  the  experi¬ 
mental  set-up  described  above.  The  first  point  to  emerge  from  these 
experiments  was  that, unlike  impacts  on  pure  PETN  where  ignition  never 
occurs  before  transparency  and  rapid  flow  of  the  sample, have  occurred, 
when  a  grit  is  present, ignition  can  take  place  before  the  explosive  has 
been  sufficiently  compacted  to  fail  plastically,  form  a  transparent  layer 
etc.  This  observation  implies  that  the  pressures  needed  to  cause  explosion 
are  much  smaller  than  with  a  pure  material. 

If  a  sample  ignites  while  it  is  still  opaque  it  means  that  small 
polymer  particles  within  the  pellet  can  not  be  seen,  and  thus  association 
of  the  initiation  sites  and  grit  particles  is  difficult.  To  overcome  this 
problem,  two  different  approaches  were  adopted.  In  the  first  instance, 
pressed  pellets  of  pure  explosive  were  surrounded  by  polymer  discs  and  the 
impact  phenomena  of  the  whole  array  photographed.  The  second  type  of  experi¬ 
ment  involved  pressing  pellets  of  explosive,  into  which  had  been  put  single 
polymer  discs,  effectively  modelling  a  single  large  polymer  particle  in  a 
sample  of  pure  explosive.  Examples  of  results  obtained  from  these  experi¬ 
ments  are  illustrated  in  Figures  8- ’ ’ • 

Figure  8  shows  a  0.65  mm  high  PETN  pellet  surrounded  by  five  I  ran 
high  PS  discs.  During  the  impact  the  PETN  disc  (dark  central  area)  breaks 
up  and  is  forced  into  close  contact  with  the  polymer.  The  polymer  discs 
deform  plastically  and  then  undergo  catastrophic  failure  as  described 
earlier.  In  frame  8b  all  the  discs  are  exhibiting  cracking  and  shear  band- 


1 1 


ing  and  by  frame  8c  have  failed  catatrophically.  In  frame  8d  we  can 
see  the  first  ignition  site  (labelled  I)  with  reaction  spreading  along 
two  paths  into  the  PETN.  A  second  site  (labelled  2)  has  developed  in 
frame  8e .  By  the  final  frame,  fast  reaction  has  spread  throughout  most 
of  the  PETN. 

An  example  from  the  second  type  of  experiment  is  shown  in  Figure  9. 

In  this  case  a  single  polymer  disc  of  PS  was  pressed  with  the  PETN 
powder  to  form  a  composite  disc.  The  selected  frame  shows  ignition  at 
three  points  (I  to  3)  on  the  polymer/explosive  interface.  The  polymer 
is  the  lighter  central  area,  and  the  darker  annulus  is  the  PETN.  In 
the  lower  part  the  PETN  is  becoming  ligher  and  this  is  due  to  the  onset 
of  fusion.  At  its  outer  edge  the  PETN  is  beginning  to  produce  fast  jets. 
The  mottled  appearance  of  the  polymer  is  partly  due  to  catastrophic  failure 
and  partly  due  to  a  very  thin  layer  of  trapped  PETN.  This  is  why  reaction 
also  spreads  inwards  as  well  as  outwards.  Frame  9b  shows  a  later  state 
where  fast  reaction  is  more  developed.  The } expansion  velocity  of  the  PS 
disc  just  prior  to  ignition  was  ca.  35  m  s  .  The  catastrophic  failure 
occurred  between  frames.  Experiments  with  PC  showed  that  it  behaved 
in  a  very  similar  manner,  with  cracking  and  rapid  radial  expansion  followed 
by  ignition  at  the  interface. 

The  situation  with  PP  is  quite  different  and,  although  much  radial 
expansion  at  quite  high  velocities  is  observed,  the  polymer  does  not  crack 
or  shear  locally  and  ignition  does  not  occur. 

The  results  obtained  with  HMX/polymer  systems  are  similar  to  those 
obtained  with  PETN.  Figure  10  illustrates  an  impact  on  a  20  mg  sample 
containing  a  PP  disc  and  Figure  11  a  20  mg  sample  containing  a  PC  disc. 

In  both  cases  the  polymer  is  the  light  central  area  and  the  explosive  the 
darker  area.  Figure  10  shows  that  yield  and  flow  of  both  the  PP  and  HMX 
occur  but  no  ignition  results.  In  Figure  11  ignition  occurs  at  the  polymer/ 
HMX  interface  (frames  d  onwards)  after  catastrophic  failure  of  the  PC. 

The  findings  of  the  photographic  work  described  in  this  section  may 
be  summarised  as:  (i)  ignition  occurs  after  or  during  the  catastrophic 
failure,  cracking  and  fast  flow  of  the  polymer;  (ii)  the  initiation  site 
is  very  closely  associated  with  the  position  of  the  polymer  and  (iii)  of 
the  polymers  tested,  only  those  which  undergo  catastrophic  failure  are 
effective  in  promoting  ignition. 

2 .-  Impact  Behaviour  in  the  Drop-Weight  Test 

(a)  Impacts  of  polymer  samples 

Samples  of  the  polymeric  materials  were  impacted  in  a  drop-weight 
machine  at  strain  rates  similar  to  those  observed  in  the  photographic 
work  (e  ca.  10  s  ')so  that  differences  in  their  mechanical  behaviour 
could  be  observed.  The  drop-weight  machine  used  was  similar  to  that  des¬ 
cribed  by  Heavens  and  Field  (1974)  and  is  illustrated  schematically  in 
Figure  12.  The  sample  was  placed  between  rollers  R1  and  R2  which  had  their 
edges  ground  down  to  give  a  sharply  defined  contact  area  of  diameter  5  mm. 
The  samples  used  were  discs  of  polymer  typically  1.5  mm  high  and  7  mm  dia¬ 
meter.  The  contact  area  was  lubricated  with  colloidal  graphite.  Two 
semiconductor  strain  gauges  were  mounted  in  series  diametrically  opposite 
to  each  other  on  R3  in  order  that  the  resistance  change  induced  in  the 
pair  of  gauges  during  impact  would  be  due  to  compressive  stresses  alone. 
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The  system  was  calibrated  statically  using  an  Instron  machine. 

As  the  sample  diameter  was  greater  than  the  contact  area,  the  area  of 
contact  during  impact  remained  constant  throughout  the  whole  process.  The 
pressure-time  traces  were  obtained  on  paper-tape  utilising  a  Data  Lab  DL922 
transient  recorder  and  analysed  on  a  HP9825  calculator  using  a  programme 
based  on  the  method  of  analysis  of  Af anas' ev  and  Bobolev  (1971).  The 
programme  numerically  integrates  the  force-time  trace  produced  by  the 
.equipment  to  yield  a  velotiy-time  trace  for  the  weight 

-1  1  r 

V(t)  -  VQ-  M  Jo  F(t)  dt  (3) 


with  VQ  the  velocity  of  the  weight  at  impact  and  M  its  mass.  A  second 
integration  gives  the  displacement  of  the  weight  as  a  function  of  time  z(t). 

z(t)  -  CfQ  V(t)  dt  (4) 

z(t)  is  the  compression  of  the  whole  system  (sample  and  anvils)  and  a  small 
amount  z  (t)  due  to  the  compression  of  the  system,  which  can  be  determined 
by  performing  impact  with  no  sample  present,  is  subtracted  from  z(t)  to 
yield  the  compression  of  the  sample. 

zg(t)  -  z(t)  -  *A(t)  (5) 

The  sample  strain  is  then  given  by 


e(t) 


In 


^  H 


2s(t) 


(6) 


with  H  the  initial  height  of  the  sample  and  e(t)  the  strain  at-  time  t. 
The  stress  at  time  t  is  given  by  division  of  the  force-time  curve  by  the 
area  of  contact  A. 


a(t)  ■  F(t)/A 

Elimination  of  t  between  equations  5  and  6  yields  the  stress-strain  curve. 

Results  obtained  with  PC  and  PP  are  illustrated  in  Figure  13  (a)  and 
(b).  The  behaviour  of  PS  was  similar  to  that  of  PC.  The  main  difference 
between  the  two  materials  is  the  evidence  of  catastrophic  failure  in  PC 
(as  shown  by  the  sharp  drop  at  e  *  1.8)  and  its  absence  in  PP.  Both  mater¬ 
ials  show  evidence  of  strain  softening  followed  by  orientation  hardening, 
but  the  strains  which  PP  is  capable  of  are  far  greater  than  those  achieved 
by  PC.  The  onset  of  orientation  hardening  also  takes  place  much  earlier  in 
PC  at  a  strain  of  ca.  0.5  as  opposed  to  the  value  of  2.0  observed  with  PP. 

At  strains  of  >2  the  thickness  of  the  PP  layer  becomes  comparable  to  the 
correction  factor  for  anvil  deflection  and  thus  the  calculated  strain 
become  increasingly  unreliable. 

The  shape  of  the  stress-strain  curve  may  therefore  be  inaccurate  at  high  strains 
These  results  confirm  the  findings  of  the  photographic  work  in  showing 
a  marked  difference  in  high  strain  rate  impact  behaviour  between  sensitising 
and  non-sensitising  materials  and  again  point  to  catastrophic  failure  as 
being  the  process  responsible  for  a  polymer's  sensitising  action. 


VTaCj»*g  i; 
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(b)  Impact  behaviour  of  explosive/polymer  samples 


The  sensitivity  of  pure  explosives  PETN  and  HMX  and  explosive/polymer 
samples  were  compared  by  determining  their  50%  height  on  the  drop-weight 
machine  outlined  above.  The  50%  height  is  the  height  from  which  the  hammer, 
if  allowed  to  fall  freely,  would  have  a  50%  probability  of  causing  initiation. 
The  method  of  determination  was  based  on  the  "up  and  down"  method  which  is 
fully  described  in  Dixon  and  Massey  (1957). 


A  large  number  of  drop-weight  tests  were  carried  out  on  explosive/grit 
samples  to  determine  (i)  which  polymers  gave  rise  to  increased  sensitiveness, 
(ii)if  the  size  of  the  grit  particles  was  important ,(iii)  the  effect  of 
variation  of  the  percentage  weight  of  grit,  and  (iv)  if  a  polymer  which  sensi¬ 
tised  one  secondary  was  effective  in  sensitising  another.  All  the  experiments 
were  carried  out  using  twenty  25  -  1  mg  samples  of  powder  which  had  been 
pressed  statically  to  10  kbar  to  give  pellets  of  5  mm  diameter  and  approxi¬ 
mately  0.8  mm  height.  In  these  experiments,  three  grades  of  grit  were  used 
which  have  been  called  coarse,  fine  and  very  fine.  In  the  'very  fine'  grade, 
all  the  particles  had  been  through  a  76  pm  mesh,  the  'fine'  grade  indicates 
that  the  grit  was  used  as  obtained  from  the  freezer-mill  with  a  maximum 
particle  mass  of  ca.  20  Pg.  The  coarse  material  was  as-received  grit  from 
AWRE  with  particles  of  mass  up  to  1  mg  and  up  to  1  mm  long.  These  reuslts 
are  set  out  in  Table  5.  Unless  otherwise  stated,  all  the  samples  contained 
20%  by  weight  of  grit 

Table  5:  50%  Heights  for  Explosives  and  Explosive/Binder  Samples 


Sample 


50%  ht  (cm) 


PETN  (pure) 


21.8  -  1.0 


PETN  +  PC  (very  fine) 

PETN  +  PC  (very  fine)  repeat 


10.5  -  2.0 
10.7  -  0.6 


PETN  +  PS  (very  fine) 

PETN  +  PS  (fine) 

PETN  +  PS  (coarse) 

PETN  +  PS  (1  mm  high,  1mm  diam.) 


11.1  -  0.3 
13.0  ±  1 .5 

10.4  ±  0.6 

10.5  ±  3.0 


PETN  +  PMMA  (coarse) 
PETN  +  PMMA  (fine) 


21.0  1  1.6 

19.2  ±  0.8 


PETN  +  PP  (very  fine) 

PETN  +  PP  (very  fine)  10%  grit 


>8.8  t 
9.2  ± 


HMX  (pure) 


HMX  +  PC  (fine) 


16.0  -  1.5 


HMX  +  PS  (very  fine) 


15.1  -  1.0 


HMX  +  PMMA  (coarse) 


22.6  -  1 .0 


HMX  +  PP  (fine) 


24.4  -  0.8 


AgN^  (pure) 


19.3  -  0.8 


6.8  -  1.3 


AgN3  +  PC  (fine) 


From  che  results  set  out,  it  can  be  seen  that  the  size  of  the  grits 
used  does  not  have  any  great  effect  on  the  results  produced.  It  might  be 
expected  that,  if  the  effectiveness  of  the  grit  depended  on  its  ability  to 
disturb  the  flow  of  che  explosive  as  i t  was  extruded  from  between  the 
anvils,  the  larger  particles  by  virtue  of  their  size  would  prove  more  effec¬ 
tive,  the  very  fine  particles  being  so  small  as  to  be  carried  along  in  the 
flow  of  the  explosive.  The  results  obtained  show  that  this  is  not  a  major 
factor . 

The  results  also  show  that  the  grits  are  not  specific  in  their  action: 

PC  and  PS  being  equally  effective  in  sensitising  HMX  and  ’’ETN,  while  PP 
and  PMMA  only  slightly  sensitise  both  of  these  secondar  ’  explosives.  As 
noted  before  (see  for  example,  Heavens,  1973,  and  Heavens  and  Field,  1974) 
primary  explosives  such  as  lead  azide  Pb(N-)0  and  silver  azide  AgN^  are  not 
particularly  sensitive  in  the  drop-weight  test  provided  the  anvils  are  well- 
aligned.  However,  as  the  table  shows,  the  polymers  do  have  a  significant 
effect  on  the  sensitivt ness  of  AgN^  and  in  a  consistent  way.  Finally,  in 
one  set  of  experiments  the  percentage  of  grit  was  reduced  to  10%  but  this 
caused  no  significant  change. 


2 . '  Thermal  Properties 

(a)  Mechanical  effects  and  thermal  properties 

The  mechanical  failure  properties  of  polymer  additives  appear  to  be 
of  major  importance  in  determining  their  sensitising  effect.  Those  mater¬ 
ials  which  undergo  bull  plastic  deformation  proving  ineffective  as  sensitisers 
and  those  which  fail  c;  tastrophically  and  with  intense  local  deformation 
effective.  Thermal  prc oerties  areof  majoi  significance  in  determining 
the  failure  modes  of  pc  Lyme rs  particularly  when  high  strain  rates  are  in¬ 
volved.  This  is  becauss  they  soften  at  relatively  low  temperatures  and  have 
low  densities,  specific  heats  and  thermal  conductivities.  Additionally, 
materials  can  desensitise  explosives  by  effectively  acting  as  heat  sinks 
and  quenching  hot  spots  before  they  can  fully  develop.  The  work  of  Bowers, 
Romans  and  Zisman  (1973)  for  example,  indicates  that  sensitivity  decreases 
with  increasing  specific  heat  of  additive,  those  additives  which  are  most 
effective  in  this  respect  being  materials  such  as "Superla  wax"or"Carbowax" 
which  have  low  melting  points  and  thus  absorb  their  heat  of  fusion. 

In  cases  where  thermal  properties  were  not  known  they  were  measured 
using  a  Perkin-Elmer  DSC-2  calorimeter. 

The  results  obtained  show  that  both  PC  and  PS  have  specific  heats 
which  range  from  1.25  kJ  kg-*  K  '  at  room  temperature  to  2.35  kJ  kg-'  K-' 
at  600  K,  the  increase  aeing  reasonably  linear.  PC  has  a  melting  point 

of  420  to  430  K  which  oily  introduces  a  very  slight  increase  iiijthe  effective 

specific  heat,  giving  a  latent  heat  of  fusion  of  ca.  7.2  kJ  kg  .  PS  shows 
melting  over  the  range  i70-500  K  with  a  heat  of  fusion  of  ca.  6.6  kJ  kg  '. 

PP  has  a  greater  specific  heat,  varying  from  1.92  kJ  kg-'  K-'  at  room  temp¬ 

erature  to  2.09  kJ  kg-'  K-'  at  450  K.  It  melted  over  the  range  390  - 
460  K  with  a  latent  hea:  of  fusion  of  ca.  64,0  kJ  kg-'  K-' .  PMMA  had  a 
relatively  constant  spe:ific  heat  of  2.9  kJ  kg-1  K-'  over  the  range  470  - 
630  K  where  a  large  dec  imposition  endotherm  was  produced.  From  the  results 
set  out  above,  it  could  be  predicted  that,  on  the  basis  of  thermal  properties 
alone,  the  order  for  ho  :-spot  quenching  for  the  polymers  investigated  would 
be  PP,  PMMA,  PC,  PS. 
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interface  to  be  viewed  by  an  infra-red  detector .  Friction  cannot  be  regarded 
as  just  a  surface  shearing  effect  since,  if  one  of  the  materials  is  harder 
than  the  other,  the  asperities  on  the  harder  surface  can  plough  out  the  softer 
material  to  an  appreciable  depth  below  the  surface.  it  is  thus  possible  for 
considerably  and  varied  deformations  to  take  place  in  a  frictional  situation, 
and  estimates  of  the  surface  temperatures  achieved  are  of  interest  both  to 
the  present  investigation  and  the  subject  of  polymer  friction  and  wear  in 
general . 

(a)  Apparatus  and  technique 

The  apparatus  constructed  to  perform  these  experiments  was  basically 
a  modification  of  that  described  by  Bowden  and  Thomas  (1954)  and  is  shown 
schematically  in  Figure  14.  It  consisted  of  a  50  mm  diameter  rotating  sapphire 
disc  (infra-red  transmitting  to  ca.  5.3  pin)  against  which  a  polymer  pin  was 
loaded  by  means  of  a  balance  arm.  The  area  of  the  interface  under  observation 
was  defined  by  a  0.8  mm  slit  mounted  directly  above  the.  rotat  ing  disc.  The 
disc  rotated  at  a  frequency  of  ca .  56  Hz  giving  an  interfacial  velocity  of 
ca.  7ms"*.  Above  this  disc  was  a  l'MMA  chopper  with  radial  regions  i.n  the 
sequence,  hi  aeked-out ,  slit,  l’MMA ,  slit,  blacked  out  etc.  The  di.se  was  ro¬ 
tated  so  that  the  radiation  from  the  interface  was  chopped  at  a  frequency 
of  10  V.llz.  This  allowed  temperature  variations  of  duration  1  ms  or  greater 
to  be  followed  (since,  for  any  one  temperature,  determination,  a  complete 
blaeked-out  to  blacked-out  sequence  must  be-  obtained,  i.e.  5  'chops’).  After 
chopping,  the  radiation  was  detected  by  a  Mullard  Rl’Y  36  infra-red  detector 
which  is  sensitive  over  the  range  1  -  5.5  pm  with  a  sharp  cut-off  at  5.5  pm. 
The  detector  was  biased  to  give  the  best  signal -to-noi se  ratio.  The  detector 
output  was  taken  to  a  Tektronix  7A22  amplifier  and  the  amplifier  output  was 
fed  to  a  Datalab  DL922  transient  recorder.  Traces  wore  then  punched  onto 
paper- t ape  for  analysis  on  the  Hewlett-Packard  HP  9825  desk-top  calculator 
and  plotter. 


Ihc  purpose  of  the  l’MMA  filters  was  to  provide  an  absorption  band 
(ca.  3.5  ym.)  which  would  absorb  a  differing  proportion  of  the  signal  depend¬ 
ing  on  the  temperature  of  tie  source.  Thus,  by  calculating  the  ratios  I/T„ 
(with  I  the  intensity  through  the  slit  and  7,,  through  the  filter)  for  sources 
of  various  known  temperatures,  a  calibration  curve  of  1/l.j.  against  temperature 
could  be  drawn  and  estimates  of  temperatures  of  other  sources  made  using  the 
curve.  In  order  for  the  calibration  to  ha  valid,  it  was  assumed  that  the 
sources  (calibration  and  specimen)  were  hi  having- as  black-bodies.  Problems 
associated  with  different  substances  having  differing  black-body  cmi ssivi t ios 
are  eliminated  by  the  use  of  the  ratio  I/lj,,. 


(b )  Cal i brat  ion 


The  theory  of  calibration  follows  the  treatment  of  Land  (1944)  and 
Parker  and  Marshall  (1943).  From  Planck's  radiation  law  we  have 

Ea  t  -  C  ex  X'5  /  <*n/AT  -  |)  (7) 

with  >y  the  energy  emitted  by  the  black-body  at  wavelength  \  and  temperature 
T,  C  and  B  constants  and  the  emissivity  of  the  surface  at  wavelength  X. 
This  can  be  approximated  to  Wien's  formula. 

K  =  A  G  c\  t  X  ^  exp(-B/A  T)  (8) 

P  P  P 

with  E  the  radiant  energy  received  by  the  detector, A  a  constant  dependent 
on  the  area  of  the  source,  C,  c  constant  dependent  on  the  geometry  of  the 
system,  t  the  transmission  coefficient  of  the  sapphire  and  X  the  effective 
wavelength  of  response  of  the  detector.  The  relation  betweeft  the  response 
of  the  detector  H  and  the  energy  falling  on  it  is  of  the  form 

H  =  Em  (9) 

with  m  a  constant. 


Thus  from  equations  8  and  9: 
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In  H  -  ro(ln  (A  C  t)  ♦  ln(c,  X  )  -  B/X  T) 

*  P  P 


00) 


When  the  radiation  passes  through  the  PMMA  filter  the  optical  properties 
of  the  system  are  changed  (I,F.,A  and  XpC.bange)  but  the  factors  A,  G  and 
t  remain  the  same.  Thus  if  H*  is  the  reduced  output  on  passing  through 
•the  filter  we  may  write 


In  H'  *  m(in(AGt)  +  in(c| 
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From  equations  10  and  11 
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This  relation  is  independent  of  the  source  area  and  geometrical  fractors  and  is 
independent  of  the  absolute  values  of  the  emissivities.  It  depends  only  on 
the  differences  in  optical  properties  of  X  and  X'  (and  this  is  dependent  only 
on  the  properties  of  the  system  with  and  w?tbout  Phe  filter).  Thus  a  plot 
of  lnjj,  against  J.  should  give  a  straight  line  and,  in  fact,  this  is  what 
is  found.  T 


Initially  calibration  vras  carried  out  using  a  2.5  mm  thick  PMMA  chopper 
disc  and  a  0.27  mm  diameter  platinum  wire  as  an  artificial  hot-spot.  The 
temperature  of  the  wire  was  varied  by  changing  the  DC  current  flowing  through 
it.  The  voltage  arrows  a  central  portion  oi  the  wire  was  monitored  via  two 
25.4  pm  diameter  platinum  wires  which  were  spot-welded  to  the  larger  wire  43.0  mr. 
apart.  These  wires  were,  thin  enough  not  to  affect  t  l.e  temperature  of  the  larger 
wire  which  was  considered  to  be  uniform  over  the  46.0  mm  long  middle,  portion. 
From  the  knowledge  of  the  resistance  of  this  central  portion,  its  temperature 
was  determined  using  standard  tables.  The  central  portion  was 


viewed  by  the*  detector  though  a  slit  of  width  1  inm  mounted  perpend  i  cular  to 
Clio  length  of  the*  wire.  Alliums',!)  this  yaw  a  gopd  1  i near  relationship  over 
the  range  770  - 1 270  C  (Figure  1 5a) ,  the  ratios  -#  for  lower  temperatures 
became  to  high  to  he  measured  accurately  and  errors  associated  wiLh 

tlie  determination  of  current  ami  voltage  across  the  wire  became  significant 
so  a  different  approach  was  adopted  for  the  low  temperature  calibration 
of  the  system. 

The  thickness^of  the  PMMA  chopper  was  reduced  to  I  mm  in  an  attempt  to 
improve  the  ratio  ^  ,  and  a  variable  temperature  soldering  iron,  to  the  tip  of 
which  was  attached  a  thermocouple,  was  used  as  the  source.  The  calibration 
produced  is  illustrated  in  Figure  15(b).  A  good  linear  relationship  is  again 
obtained  until  the  temperature  drops  below  ca.  500  If  (-,  =  2.1  x  10~^)  when 

the  seat  ter  in  the  results  begins  to  get  quite  considerable.  This  is  because 
the  ratio  has  again  increased  to  values  where  the  accuracy  in  its  determination 
is  reduced  by  noise.  At  such  low  temperatures,  the  power  emitted  by  a  source 
which  is  in  the  range  detectable  by  the  detector  is  very  low  and  eonsiderat ions 
of  si rnal -to-noise  ratio  is  to  increase  the  size  of  the  source.  A  simple  Method  . 
of  improving,  the  signal  to  noise  ratio  is  to  increase  the  size  of  the  source . 

This  has  been  tried  but  there  are  problems  associated  with  doing,  so  which  will  be 
described  below.  Further  reduction  in  the  thickness  of  the  I’MIA  filter 
would  improve  the  rat  io  -  i  >  however  1  mm  is  the.  smallest  commercially 
available  thickness  and  it  would  be  necessary  to  use  other  materials  to 
improve  the  system  in  this  manner.  In  fact,  0.2  mm  thick  acetate  sheet  was 
tried  but  was  found  to  be  unsuitable  as  the  position  of  its  absorption  bands 
was  sue!)  that  very  little  change  in  the  value  of  II'  was  produced  for  a  wide 
range  (400-600  K)  in  T. 

The  calibration  v.ns  checked  by  measuring  the  radiation  received  from  a 
t hormone l or  bulb  placed  beneath  a  sapphire  disc.  The  bulb  was  heat  cd  by  a 
soldering  iron  and,  with  the  iron  removed,  a  temperature  reading  was  taken, 
and  the  result  obtained  from  the  infra-red  system  and  the  thermometer  scale 
compared.  For  a  thermometer  reading  of  4 8f>  K,the  infra-red  system  gave  a 
value  of  489  K. For  453  K  the  system  gave  438  K  and  for  the  low  temperature 
of  383  K  the  system  gave  418  K.  It  can  thus  be  seen ,  and  it  is  clear  from 
the  calibration  curve,  that  the  determination  of  temperatures  below  ca.  500  K 
leads  to  an  increasing  inaccuracy  in  the  result.  In  the  results  quoted  1  at  ti¬ 
the  error  in  the  result  is  estimated  by  drawing  an  envelope  which  encloses 
all  the  calibration  points  for  values  of  ,J,  less  than  2.1  x  10  (dotted 
in  Figure:  15b).  The  value  of  489  K  for  an  'actual'  temperature  of  488  K  may 
be  rather  fortuitously  close  and  all  estimates  of  temperature  will  have  an 
error  of  the  order  of  *  10  k. 

The  preliminary  analysis  of  the  results  was  carried  out  by  the  HP  9825 
calculator.  A  programme  to  perform  this  analysis  was  written  which  changed  the 
zero  level  of  the  traces  such  that  it  intersected  the  pcak-to-pcak  levels 
approximately  half-way  (see  Figure  16).  The  heights  of  the  peaks  were  then 
dedu« cd  by  finding  the  ha  1 i -way  points  be  tween  the  intercepts  of  the  trace  and 
the  axis,  taking  as  the  numerical  value  of  the  peak  height  an  average  of  seven 
point  s  around  (an  !  im  lading)  this  hal f-vuy  value.  The  trace  was  then  pi  "l  t  cd 
out  teg,,  thcr  will,  a  list  in,;  of  the  pea1;  heights  ana  the  remainder  of  the  ana¬ 
lysis  (finding  the  ratio  tj,,  v  is  done  bv  b.md\  Figure  16  shows  three  typical 
result,:,.  ( i )  is  for  a  soldering  iron  at  '->75  k  ( i  i )  for  a  small  diameter  PC  pin 

am'  (ii.i)for  a  large  diameter  I'C  pin.  The  improvement  in  s  i  gnal  tr  -noi  se  )  ."  t  in 

with  a  larger  diameter  pin  is  clear  from  n  cempca  i  son  of  ( i  i )  and  (tij)  In  tb 

of  v.nieh  give  approximate  I  y  tin-  same  mini  inn  value  ot  the  ratio 
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(c)  Results 


i<>. 


In  the  experimental  investigation,  polymer  pins  of  two  different  dia¬ 
meters  were  used:  small  pins  of  diametei  ::a.  1.5  nun  and  larger  pins  of  dia¬ 
meter  ca.  6  nun.  The  smaller  diameter  pins  were  used  in  the  initial  experiments 
but  for  the  later  measurements  large  pins  were  used  since  the  smaller  pins 
gave  a  poor  signal-to-signal  noise  ratio  and  wore  very  easily  deformed  due 
to  their  large  length/diameter  ratio.  The  main  problem  with  the  use  of  the 
larger  pins  was  that  considerable  heating  of  thesipphite  disc  took  place 
and  surface  temperatures  of  up  to  360  K  were  recorded  on  placing  a  thermo¬ 
couple  in  contact  witn  the  disc  x mud  lately  (1-2  sj  after  doing  a  run. 

Due  to  the  considerable  r derma l  mass  of  the  disc,  its  rate  cf  cooling  was 
quite  slow  (!C  -  20  degrees  per  minute!  and  sa  the  figure  quoted  above  gives 
a  good  indication  of  the  bulk,  temperature.-  reached  by  the  sapphire  in  the 
course  of  an  experiment. 

An  experiment  in  which  the  signal  level  was  recorded  with  the  large  pin 
in  contact  with  the  disc  »td  again  : -vm'dia < ■  1  y  niter  its  removal  indicated  that 
the  energy  due  to  the  bulk  heating  of  the  di«c  accounted  for  approximately 
one  third  of  the  tore'  signal.  Such  a  1  a>-pe  contribution  from  a  lower 

temperature  ,->our<e  will  give  a  mixed  signal  which  will  tend  to  reduce  the 
measured  black -body  temperature  to  a  value  near  the  mean.  Thus  the  readings 
obtained  with  the  larger  pins  will  probably  be  underestimates  of  the  interface 
temperature.  In  the  table  of  results  given  l”- low  (Table  7)  the  maximum  black- 
body  temperature  recorded  by  the  system  for  any  given  polymer,  irrespective  of 
the  typr  of  pin  used  is  quoted. 

The  results  will,  of  c.-urso  .  be  averages  over  the  total  portion  of  the 
interface  area  vithi.  the  detector’s  fi«.*ld  of  view  and  it  is  quite  likely  that 
some  regions  (those  being  highly  deformed)  will  have  temperatures  above  this 
average  value  and  otlu  rs  will  have  lover  values.  by  far  the  largest  amount 
of  experimental  time  vu  ■.  sport  -n  the  polymers  FV  .and  PC  re  that  with  an 
increased  amount  of  ohn rvat ion  time  it  may  be  possible  to  record  higher 
maximum  value s  of  tempera  faro  for  the  « hree  other  material s .  it'  is,  however,  in¬ 
teresting  that  the  reraj  •.  rnlvrcj  obtained  follow  the  order  obtained  for  the 
sensitiveness  of  tin.  pel  vsner/exp!  a  •  vo  mixtures.  It  is  also  clear,  particularly 
in  the  case  of  VC,  Mint  temperature  we  1 1  above  the  commercially  quoted 
softening  temperature  rtf  4  10  K  and  tin-  observed  ’melting’  temperature  (in 
previous  work)  of  4  30  K  nr."  obtainable.  In  the  case  of  PP  however  this  wan  nor 
found  to  be  the  case  and  no  tempers! ures  above  the  commencement  of  melting 
at  390  K  wore  obtained. 

Tabic  7:  Tempo Tarures Ac  In  eyed  During  Polymer  Friction 


Pol  ymer  Maximum  Intt-rfacial  Tcmpcrat ure/K 


4 

pp 

360  -  30 

I’MMA 

420  -  30 

«fr 

PS 

430  -  20 

■4* 

F’C 

490  -  10 

20. 


I* 

* 

L  •*  ' 
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Kxainln.-'t  ion  «>l  the  polymer  pins  after  flic  frict  i.'nal  exper  imem  s  reveal 
that  in  flu-  ease  of  10  and  I’S  severe  do format  ion  of  the  polyim-r  bad  been 
prediu  ,•(!  with  evidence  of  much  si  ret  chi  nr.  and  shearing  of  molecular  layers. 

The  1*1*  pin,  however,  present  a  much  more  rounded  appearance  which  indicated 
that  melting  of  the  surface  layers  had  taken  place. 

Figure  17  shows  the  overall  appearance  of  the  pins  after  the  experiments 
(a)  is  a  PC  pin  and  the  many  layers  which  have  sheared  off  and  built  up  on 
the  trailing  edge  are  clearly  visible,  (b)  is  a  PP  pin  and  shows  the  rounded , 
melted- 1  coking  appearance  of  the  surface,  (c)  is  a  PMMA  pin  and  shows  the 
roughened  surface  obtained  with  no  build  up  on  the  trailing  edge.  All  the 
pins  in  Figure  17  were  of  diameter  6  mm. 

Microscopic  examination  confirmed  that  very  severe  deformation  of  PC  had 
taken  place  with  clear  evidence  of  considerable  strain  in  the  frictional 
direction  (.Figure  18(a)).  PMIaA  (Figure  18(h))  presented  a  rather  *  powdery ' 
appearance,  with  many  particles  of  material  on  the  surface.  The  wearing  away 
of  the  P.-.MA  pins  proceeded  in  a  more  hri  t  tie  manner  wi  tli  small  particles  chipped 
out  of  Lhr  surface  and  with  no  accumulation  of  sheared  or  melted  layers  on  the 
trailing  edge  of  the  pin.  The  .surface  of  PP  on  the  other  hand  was  relatively 
featureless;  as  might  bo  expected  from  a  me!  tied  and  refrozen  layer. 

Although  there  is  clear  evidence  that  PC  hid  deformed  by  the  shearing  off 
of  layers,  there  was  no  evidence'  of  catastrophic  failure  having  taken  place: 
the  shearing  could  be  followed  by  the  naked  eye  and  took  place  over  a  time- 
scale  ot  cn.  1  s  for  a  6  mm  pin.  Nonetheless,  temperatures  well  above  the 
consnturioal  softening  temperature  were  recorded  with  PC,  and,  vlu.n  allowance 
is  made  Cm  t lie  fact  that  these  were  averages  over  the  whole  surface,  it  does 
seem  pi"-.:' i  hi  c  that,  under  conditions  of  extreme  (and  probably  locally  rapid) 
deformation,  very  high  temperatures  could  be  achieved  by  polymers  such  as  VC. 
These  temperatures  would  he  well  above  the  limits  set  by  the  old  ideas  of  the 
temperature  rise,  achieved  in  friction  being  limited  by  t  lie  melting  point  of 
the  lower  molting  point  material  which  was  based  on  the  frictional  properties 
of  metals.  Any  increase  in  melting  point  of  the  polymers  due  to  the  force 
applied  by  the  counterweight  would  be  negligible  since  the  maximum  weight 
used  with  the  6  mm  pins  was  ?.  kg  and  wi  th  the  1.5  mm  pins  was  500  g  giving  a 
maximum  stress  of  3  IJI’a,  which  corresponds  to  ca.  30  bar.  This  would  have  an 
insignificant  effect  on  the  melting  point,  variations  of  ca.  200  bar  being 
needed  to  cause  an  appreciable  (ca.  10  K)  change  (Zollor  I97S).  It  does 
therefore  appear  that  high  temperatures  can  he  reached  during  the  deformation 
of  polymers,  and,  in  the  design  of  any  system  containing  a  polymer,  considernti 
must,  be  given  to  the  possibility  of  high  temperature  regions  during  deformat  ion 

2.7  Com  1  us  i('_ns  and  )hi  sens  s  i  on 

A  combination  el  high-speed  photographic,  pressure  measuring  and 
thermal  analysis  t hit  i  cues  have  been  used  to  study  the  behaviour  and  inter¬ 
action  et  i-::p]nsi  v-  ai-d  poly,-.,  rs .  Chemical  interactions,  for  the  materials 
tested,  were  found  to  be  of  min-';  imp.  rt.uce  compared  with  mechanical  cons  id¬ 
em  l  1 1  n  •  is  imp  a,  t  sea:  itis-iti,  ■  .  Pure  explosives  wore  found  to  ignite  only 
after  ivp  ,1  .inn  m  •  re  do  f  <  ri’in :  inti  of  the  s.-.trpl  ••  bad  occurred.  In  the 
pres,  i:  ,-f  ?  or,  .  i  :  i  is;  -,r  it  s  >■ ,  n-r  ,  Ml,  severe  d.-rei  ion  was  not 

0  e  ■  - 1  •  •  y  !,r  s  or. .  hot'1  :  r>  ii.'.etur,  and  the  rim*  I  proper  t  i  of  >’f  an 

-  '.hi  ’  •  i  is  j  :  Y  ■  :  .  i  i  in  f  i,"  .lie.;  i  !  •:  s  •  -.  ,  >  s  i  l  i  \  l  v  .-  f  i  ,.e  t  .  A  p.  ,1  yrx  r 
b  ■'  i  rut  : ,,f  hi  ;  !i  with  1  ■  >w  spec  i  f  i  e  h  .-at  and  lev  li--.it 

'I  Ins  r,  -'b  ■  e  t  ,  i  •  i  -  i-,,|  ii.-i  i  ,  Tui'  f  a  i  !  n  I  1  be  likeiy  to 


cause  sensitisation.  In  order  to  tost  this  hypothesis,  two  polymers  were 
chosen  on  the  basis  of  their  strength  and  thermal  capacity;  polyester  which 
has  high  strength  and  low  capacity  and  high  density  polyther^  (HOPE)  which 
has  a  low  strength  and  high  capacity.  They  were  mixed  with  PETN  and  tested 
in  the  drop-weight  machine  in  the  same  manner  as  the  samples  described  earl  ier. 
The  results  obtained,  together  with  a  few  comparison  results  from  Table  5  are 
set  out  in  Table  8.  Thermal  properties  can  be  found  in  Table  6. 

Table  8:  50%  Weights  for  Expl osive/Polymer  Samples 


PETN  (pure 


PETN  ♦  HDPE 


PETN  +  PP 


PETN  +  Polyester 


PETN  +  PS 


50%  height  (cm) 


21.8  -  1.0 


24.8  -  ) .5 


18.8  -  0.8 


15.5  -  1.0 


11.1  -  0.3 


These  results  confirm  the  predictions  that  polyester  should  sensitise 

and  HDPE  should  not.  In  fact  it  appears  that  HDPE  is  acting  as  a  desensitiscr . 

The  conclusion  reached,  therefore,  is  that  polymers  which  are  most 
likely  to  sensitise  are  those  such  as  PC  which  are  tough  below  their  glass 
transition  temperature  (as  opposed  to  PMMA  which  behaves  in  a  much  more 
brittle  manner),  and  are  likely  to  fail  catastrophically.  Those  least 
likely  to  sensitise  are  those  which  have  a  low  glass  transition  temperature 
(PP  or  HDPE)  and  which  deform  in  bulk  with  little  local  concentration  of 
energy.  Polyester  is  somewhat  intermediate  ir.  behaviour  (and  sensitising 
effect)  between  PC  and  PMMA  as  it  it  tougher  than  PMMA  but  much  more  brittle 
than  PC,  and  so  fits  quite  well  into  the  expected  pattern. 


In  the  ease  of  the 
important  property  (Kowd 
tially  because  they  are 
relatively  large  latent 
point.  In  contrast  the 
a  range  of  temper. at  ores 
If  heat  is  produced  hv  loo 
However,  when  softening 
maintain  a  high  visrof.it 
tribute  further  temporal 
losses  in  the  case  of  a 
materials  because  >  f  the 
and  the  hot  spot  is  not 


former  the  melting  point  of  the  particle  is  the 
on  and  Curt  on  1949).  That  this  is  so  is  cssen- 
materials  with  welL-defined  melting  points  and 
heats:  the  hot  spot  is  generated  to  the  melting 

polymeric  materials  which  sensitise,  soften  over 
and  have  low  latent  hears  and  thermal  conductivities, 
a lised  catastrophic  failure  hot  spots  will  form, 
of  the  polymer  starts  the  material  will  still 
y.  Continuing  deformation  will  be  able  to  r.on- 
ure  rises  by  viscous  heating.  Thor  thermal 
polymer  are  very  low  compared  with  roost  other 
low  latent  hear,  thermal  capacity  and  conductivity, 
quenched  so  readily. 


The  friction  experiment  gave  evidence  of  temperatures  well  above  the 
softening  point.  In  this  experiment,  an  average  temperature  for  a  contact 
area  of  ca.  25  mr:~  was  maintained.  Local  hot  spot  temperature  would  be 
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expected  to  exceed  the  values  recorded.  In  other  experiments  in  the  laboratory, 
the  temperature  rise  at  the  tip  of  a  crack  propagating  in  a  polymer  was  re¬ 
corded,  again  using  I  R  techniques  (Fuller,  Fox  and  Field  1973).  Two  impor¬ 
tant  results  emerged  from  this  work.  ^T1k  first  was  that  temperatures  well 
above  the  softening  points  were  observed.  Secondly, that  the  fracture  surface 
energy  increased  with  c.ruck  vcLociLy.  This  latter  result  was 

explained  by  the  high  temperatures  causing  softening  of  the  material  at  the 
crack  tip  thus  requiring  more  energy  to  be  expended  in  crack  growth. 


It  is  clear  from  the  earlier  work  and  the  present  work  with  explosive/ 
polymer  samples  that  some  polymers  can  reach  very  high  local  temperatures 
during  their  deformation.  This  is  clearly  relevant  to  the  safe  handling 
of  explosives  since  polymeric  materials  may  be  used  for  packaging  or  as 
binding  materials  in  explosive  composites.  The  results  also  have  interest 
to  studies  of  the  fracture,  friction,  erosion  and  wear  of  polymers. 
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FIGURE  CA PTI  ON S 

1.  Activation  energy  plot  for  the  thermal  decomposition  of  PETN  for 
various  values  of  reaction  order  n.  'I he  n  value  which  gives  the  best 
straight  lino  segments  is  one  appropriate  to  the  reaction  (here  0.8). 

These  two  steps  (see  clashed  line)  corresponding  to  sublimation,  stage  I, 
and  decomposition  stage  11. 

2.  Experimental  arrangement  at  instant  of  impact,  W;  drop-weight;  G, 
glass  blocks;  M,  mirrors;  S,  sample.  The  upper  glass  block  is 
attached  to  the  weight. 

3.  Impact  on  a  disc  of  PP.  The  sample  deforms  by  bulk  plastic  flow 
throughout.  Diameter  of  field  ot  view  20  mm.  Frame  tirne/ps  : 

a, ();  l>,70;c,26b.  The  zero  times  given  in  this  and  subsequent  sequences 
do  not  correspond  to  time  after  initial  contact. 

4.  Impact  of  a  disc  of  PS.  The.  sample  deforms  initially  by  bulk  plastic 
flow  (frames  a  and  b)  but.  then  deforms  catastrophically  (frame  c) . 

Diameter  of  field  of  view,  2.5  mm.  Frame  times/ps  :  a,0;b,2l  ;c,28. 

5.  Diameter  versus  time  traces  for  impacts  on  PP  and  PS.  Mote  that 

the  diameter  of  the  1*S  has  approximately  doubled  before  the  catastropic 
failure,  but  that  the  PP  disc  has  a  final  diameter  several  times  its 
initial  value. 

6.  Micrographs  of  a  deformed  PS  disc  showing  that  the  catastrophic  failure 
is  a  result  of  fracture  and  localised  shear.  The  enlarged  view  in  (b) 
is  from  region  S  in  (a).  The  view  has  been  rotated  no  that  the  shear 
bands  arc  horizontal. 

7.  Impact  on  a  pellet  of  PETN.  Frame  a  shows  the  5  mm  diameter  pellet 
before  impact.  Rapid  flow  and  jetting  occur  as  for  powdered  layers. 

Fusion  of  the  layer  starts  in  frame  d.  The  whole  layer  is  transparent 
and  flowing  rapidly  by  frame  g  when  the  first  to  two  ignition  sites 
forms.  The  circular  feature  in  frame  h  is  an  artifact.  Diameter  of 
field  of  view  20  mm.  Frame  times/ps  :  a ,0 ; b ,67 ;c , 22 1 ;d , 308 ;e , 389 ; 

f  ,469;g,476;h,482;i4p.9. 

8.  PETN  pellet  surrounded  by  5  PS  discs  Set:  text  for  details.  Two  ig¬ 
nition  sites  at  the  polvmer/explosive  interface  are  labelled  in  frames 
d  and  e.  Diameter  of  field  of  view  20  mm.  Frame  limes/ys  :a,0; 

b ,  9 1 ; c 1 40 ; c , Iff ; c , I t>8 ; f ,  i 75 . 

9.  25  mg  disc  of  PETN  containing  a  PS  disc.  Three  ignition  sites  are 

labelled.  Field  of  vuew  20  mm.  Frame  tiirfcs/ys:  a,0;b,14. 

10.  A  PP  disc  (light)  pressed  into  an  HMX  pellet.  The  first  frame  is  before 
impact.  Considerably  bulk  deformation  of  the  PP  occurs  but  there  is 

no  ignition  of  the  explosive.  Field  of  view  20  mm.  Frame  times/ps: 
a , 0 ; b , 9 1 ; c  ,  141  ;d ,  2  14  ;o ,247 ;  f  ,402  . 

11.  Similar  to  figure  10  but  with  a  PC  disc:  except  that  impact  has  already 

occurred  in  frame  i.  The  PC  fails  catastrophically.  The  first  of  several 

ignition  sites  can  be  detected  in  frame  d.  Field  of  view  20  win.  Frame 

t  ir.it  / ti :.  ;u,U  ;b ,  34  ;c  ,94  ;d  ,  I I  ;t  ,  1  35  ;  f  ,  147 . 


12.  Experimental  arrangement  tor  obtaining  pressure-time  curves.  W,  drop- 
weight;  H,  Ri,R2,R3,  hard  steel  rollers;  C,  cylindrical  guiding  sleeve 
S,  sample. 

13.  Stress-strain  curves  for  samples  of  (a)  PP,  (b)  PC.  Catastrophic 
failure  of  the  PC  sample  occurs  at  a  strain  of  ca.  1.8.  No  discontin¬ 
uities  are  observed  in  the  PP  curve. 

14.  Experimental  arrangement  for  friction  experiments.  D,  detector; 

S,  slits;  CD,  chopper  disc;  RD,  sapphire  disc;  PP  polymer  pin; 

M,  balance  arm;  W,  counterweight;  DS,  drive  shaft. 

15.  Calibration  curves  for  the  friction  apparatus  (a)  high  temperature 
region  (b)  low  temperature  region. 

16.  Examples  of  the  traces  obtained  with  the  apparatus  (a)  is  a  soldering 
iron  at  575  g  (b)  a  small  diameter  PC  pin  and  (c)  a  large  diameter 
PC  pin.  H  is  the  sigrral  strength  when  there  is  no  filter  between 
source  and  detector  and  H'  the  strength  when  a  PMMA  filter  is  inter¬ 
posed  between  source  and  detector. 

17.  Polymer  pins  after  friction  experiments.  (a)  a  PC  pin  showing  many 
layers  which  have  been  sheared  off  and  built  up  on  the  trailing  edge 
(b)  a  PP  pin  which  has  rounded  and  shows  evidence  of  melting  (c)  a 
PMMA  pin  which  has  a  roughened  surface  caused  by  small  particles 
being  brittly  chipped  out  of  the  surface. 

18.  Micrographs  of  polymer  pin  surfaces  after  friction  experiments. 

(a)  PC  and  shows  considerably  strain  in  the  frictional  direction. 

(b)  PMMA  which  has  a  powdery  appearance  with  many  small  particles 
chipped  out  of  rhe  surface. 
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The  remaining  three  sections  of  this  renort  (3,4,5)  are  on  reaction 
kinetics  and  their  analysis.  All  have  been  published  in  the  open  literature. 

For  completeness,  reference  i3  given  to  a  paper  (Chaudhri  et  al,  1977) 
based  on  work  performed  by  Dr.  Chaudhri  when  he  visited  Picatinny,  and  to 
other  papers  published  by  the  gr  ;nce  ,r,7'.'. 

M.M.  Chaudhri  and  J.E.  Field,  "Fast  Decomposition  in  the  Inorganic  Azides", 
from  Energetic  Materials,  Vol.  !,  Pd.  H.D.  Fair  and  R.F.  Walker  (Plenum, 

1977). 

J.T.  Hagan  and  M.M.  Chaudhri,  "Fracture  Surface  Energy  of  High  Explosives 
PETN  and  RDX",  J.  Mater.  Sci.  n  (1977)  Letters  1055-58. 

M.M.  Chaudhri,  W.L.  Garrett,  0.  Sandus  and  N.  Slagg,  "The  High  Velocity 
Detonation  of  Single  Crystals  of  Alpha-Lead  Azide",  Propellants  and  Explosives, 
2,  91-93  (1977). 

R.  Patel  and  M.M.  Chaudhri,  "Differential  Scanning  Calorimetric  Studies 
of  Decomposition  of  Beta-Lead  Azide",  Thermo.  Chinrca.  Acta.  25,  247-251, 
(1978). 

T.B.  Tang  and  M.M.  Chaudhri,  "Dielectric  Breakdown  of  Ionic  Azides",  Conf. 
on  Elect.  Insulation  and  Dielectric  Breakdown,  Washington,  (1979). 

T.B.  Tang  and  M.M.  Chaudhri,  "Dielectric  Breakdown  by  Electrically  Induced 
Decomposition",  Nature  282,  5734,  54-55  (1979;. 

T.B.  Tang  and  M.M.  Chaudhri,  "The  Thermal  Decomposition  of  Silver  Azide", 

Proc.  Roy.  Soc.  A369,  83-104,  (1979). 

T.B.  Tang,  "Analysis  of  Isothermal  Kinetic  Data  f rom  Mu] ci-Srage  Decomposition 
of  Solids?',  Thermochimica  Acta,  41,  133-136  (1980). 
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AllNIKA'l 

The  paper  describes  a  graphical  computer  method  for  analysing  TG  and  DSC 
traces  which  gives  all  three  reaction  parameters  (I),  A ,  and  /;)  characterising  an  nth- 
order  reaction  from  a  single  trace.  If  the  reaction  proceeds  in  multiple  stages  (A,  A 
and  »)  triplets  can  be  obtained  for  each  stage.  The  computer  programmes  are  basically 
simple  and  use  little  computing  time  (typically  a  few  seconds).  The  advantages  of  the 
approach  over  earlier  methods  ate  discussed.  As  a  test  of  the  method,  results  on  the 
dehydration  of  calcium  oxalate  monohydrate  are  described. 

INTRODUCTION 

The  deduction  of  the  re.w.ton  parameters  £  (activation  energy),  A  (frequency 
factor),  and  n  (reaction  order)  from  TG.  and  DSC  traces  is  clearly  of  interest.  A 
number  of  methods  have  been  used  in  the  past  to  do  this  but  all  have  unsatisfactory 
features.  In  the  work  described  here  a  number  of  Fortran  programmes  were  written 
to  extract  all  three  reaction  parameters  of  an  nth-order  reaction  from  a  single  TG,  or 
DSC  trace.  The  programmes  were  run  on  the  Cambridge  IBM  370  165  computer, 
and  the  plotter  was  used  for  graphical  output.  So  that  the  advantages  of  the  present 
approach  can  be  appreciated  earlier  methods  are  briefly  discussed. 

t-ARLIIR  mi  rttoos 

Murray  and  White 1  have  shown  that  the  temperature  at  the  reaction  peak  of 
DTA  traces  can  be  used  to  obtain  £  and  .f  from  a  number  of  traces  if  a  first  order 
reaction  (eqn.  I )  is  assumed,  where  w  is  the  fractional  residual  weight,  /  the  time,  T 
the  temperature  and  R  the  gas  constant. 


If  eqn.  (1)  is  differentiated  to  obtain  the  maximum  for  du/dr  then  cqn.  (2)  can 
be  derived,  where  $  =  d7,dr  is  the  heating  rate  and  Tm  the  temperature  at  the  rate 
maximum. 


In 


£  J_ 
R  Tm 


(2) 


A  plot  of  I /Tm  versus  In  (<t>iT*t)  yields  E'R  as  the  slope  and  In  (ART)  as  the  intercept 
of  the  straight  line  with  the  ordinate.  The  main  disadvantage  of  this  method  is  that 
it  only  uses  a  single  point  (£„)  from  a  trace  to  produce  one  point  on  the  activation 
energy  plot.  A  number  of  experiments  are  needed  to  obtain  a  reliable  value  for  £. 
There  is  also  the  limitation  that  the  theory  only  works  for  a  first-order  reaction. 
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Kissinger*  snowed  that  the  technique  can  be  extended  to  an  /itii-order  reaction 
by  deriving  cqn.  (3). 


=  AnK~* 


(3) 


He  shows  in  his  paper  that  the  product  n  tej^ 1  is  independent  of  the  heating 
rate  (</;)  and  almost  equal  to  unity.  Therefore,  cqn.  (3)  effectively  reduces  to  cqn.  12). 
The  approach  again  makes  poor  use  of  the  information  contained  in  a  DTA  trace. 
Ottawa's  analysis3  of  TG  traces  begins  with  a  general  form  of  reaction  equation. 


div 

dr 


(4) 


This  is  integrated  by  Ozawa  who  shows  that  the  solution  to  the  right-hand  side  of  the 
equation  is  a  polynomial  in  the  variable  E/RT.  Equation  (5)  is  an  approximate 
formula  which  relates  the  heating  rate  <f>  with  E. 


log  4>>  +  0.4567  -~L_  =  log  4>i  +  0.4S67  -  (5) 

where  7",  and  r2  coriespond  to  the  same  fractional  decomposition  («•)  of  the  sample. 
Since  a  plot  oflog  <f>  versus  1/7" can  be  obtained  for  each  value  of  w,  Ozawa  aggregates 
all  these  results  in  a  inaster-curve  with  a  gain  in  accuracy  of  the  final  result.  Equation 
(5)  only  assumes  that  the  function  f(tv)  is  constant  for  a  given  value  of  w. 

This  method  also  needs  a  number  of  experiments  at  different  heating  rates  to 
generate  enough  data  for  a  sin, do  activation  energy  plot.  Although  the  method  was 
developed  for  thcrmograviineiric  fata,  it  can  be  applied  to  DSC  data  by  integrating 
the  peaks  to  give  values  of  u\  A  theoretical  analysis  of  such  a  process  is  given  by  Reed 
et  a!.4. 

Bcrchardt  and  Daniels4  describe  a  method  which  allows  E,  A  and  n  to  be 
determined  from  the  shape  of  a  DTA  curve.  Their  appioach  was  applied  to  reactions 
in  the  liquid  phase.  The  type  of  reaction  to  which  this  method  is  applicable  must  have 
a  single  rate  constant  and  the  activation  energy  must  not  vary  with  temperature. 
By  considering  the  heat  transfer  equation  of  the  calorimeter  cqn.  f6)  was  obtained 


d W 
df 


*<L 

KA' 


c^I  +  kat) 
P  df  J 


(6) 


where  N  is  the  number  of  motes  present,  ,V0  the  initial  number,  K  the  heat  transfer 
coefficient,  A  the  area  under  the  DTA  curve,  Cp  the  heat  capacities  of  the  two  liquids 
and  AT  the  temperature  difference  between  the  two  cells.  Integration  of  cqn.  (6)  with 
respect  to  time  gives 


*  =  *°-  £jr  (C,  AT  4  Ka} 

where  a  is  the  area  under  the  DTA  curve  up  to  the  present  temperature. 
The  expression  for  the  rate  constant  k  of  a  reaction  of  order  n  is 

j.n  _  y<n  - 1 1  dN/df 

Nm 

where  V  is  the  volume.  Substituting  cqns.  (C.)  and  (7)  for  d,Y/df  and  W  gives 
r  ddT  r  tCAT 

k  m  f  km  7-  ’  c«* _ dr  l  kAT 

\  A'o  i  {A'(d'-n)  -c'ir}" 


(7) 


(8) 


(9) 
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In  order  to  obtain  reaction  parameters,  Borcliardt  anti  Daniels’  assumed  a  value 
for  n.  This  allowed  k  to  bo  calculated.  A  plot  of  In  k  versus  1/7  (activation  energy 
plot)  yields  a  straight  line  if,  and  only  if.  the  eoireci  value  of  n  has  been  assumed, 
thus  values  of  F,  A  and  n  can  be  obtained.  The  matn  difficulty  of  this  method  is  the 
number  of  plots  that  have  to  be  made  for  each  value  of  n  until  the  points  lie  on  a 
straight  line. 

A  number  of  computer  programmes  are  presented  by  Benin  ct  al.6  and  used  to 
solve  three  problems,  namely  (i)  the  production  of  thermograms  given  the  reaction 
parameters  E,  A,  n  and  given  that  the  reaction  is  of  the  nth-order  tyjie,  (ii )  the  calcula¬ 
tion  of  the  (E,  A,  n )  triplet  from  an  experimental  thermogram  —  this  is  the  inverse  to 
problem  (i);  and  (iii)  the  determination  of  the  reaction  parameters  and  the  complex 
reaction  mechanism  from  a  given  experimental  thermogram.  The  first  problem  is  a 
straightforward  function  plot  of  the  solution  to  the  nth-order  reaction  equation. 
The  second  problem  is  more  difficult,  and  presents  itself  to  anyone  who  wants  to 
obtain  reaction  parameters  from  experimental  data.  Benin  ct  al.6  treat  it  as  a  problem 
in  curve  fitting.  They  minimise  the  discrepancy  between  the  experimental  data  and  the 
calculated  function  by  taking  as  the  discrepancy  criterion  the  maximum  deviation  of 
the  function  from  the  data.  The  solution  of  the  third  problem  is  an  attempt  to  reveal 
the  mechanism  of  a  chemical  reaction.  The  programmes  become  quite  involved  as 
they  work  out  activation  energies  for  each  section  of  the  thermogram  by  using  a  plot 
of  reaction  rate  versus  IfT  for  a  constant  amount  of  conyeision,  a,  but  different 
heating  rates. 

The  most  general  form  of  the  fit  function  assumes  a  decomposition  process 
which  proceeds  in  a  numbci  of  individual  stages.  Each  stage  consists  of  an  elementary 
nth-order  reaction,  i.c.,cqn.  (1)  with  A,  F.  and  «•  replaced  by  A,.  and  wf, respectively, 
where  *V  denotes  the  variable  belonging  to  the  /'lli  stage  of  tlxc  reaction.  The  program¬ 
me  which  optimises  the  fit  of  this  function  to  the  data  takes  up  to4  h  of  computing  time. 
This  makes  the  method  rather  expensive.  Another  weak  point  is  the  discrepancy 
criteria  between  the  data  and  the  tit,  which  is  taken  as  the  maximum  difference.  It  is 
well  known  that  the  data  points  towards  the  end  of  a  decomposition  process  become 
less  reliable  due  to  impurities  and  zero  errors.  But  these  points  have  the  same  weight 
in  the  analysis  as  all  the  others. 

Schcmpf  ct  al.7  describe  a  programme  which  fits  a  least  square  polynomial  to 
the  weight-loss  curve.  The  quantity  dw/di"  s  calculated  by  dift'erentiation  of  the 
polynomial.  A  straight  line  is  then  filled  to  a  plot  of  log  k  versus  1/71  and  the  least 
square  criterion  is  used  again.  The  reactioi  order  is  assumed  to  be  one,  and  no 
attempt  is  made  to  eliminate  the  deviations  I  root  an  ideal  reaction  at  the  beginning 
and  end  of  the  weight-loss  curve.  Our  w  ork  sl  ow  s  that  this  is  very  critical  and  can  be 
obtained  from  an  Arrhenius  plot  performed  >y  the  computer. 

NEW  METHOD 

The  method  developed  was  based  on  the  ide:  s  of  Borchnrdt  and  Daniels'.  Our 
data  were  obtained  from  both  a  thermobalancc  (St  mton  ReJcrgft)  and  a  differential 
thermal  calorimeter  (Pei kin- Elmer  DSC-2). 

The  analysis  assumes  a  reaction  which  -  m  be  described  by 


Since  the  heating  rate  h  (deg  min'  1  is  constant,  the  time  variable  can  be 
eliminated  fiont  cqn.  (10)  by  putting  d7  "dr  -  h/60  (deg  sec' '). 
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Taking  natural  logarithms  on  both  sides  gives 


.  ,,  .  ,  A,  x  60 

In  l,<ic)  -  In  — y- 


Nofc  that  (  ~  (Jii’/dr)  is  ;i  positive  quantity  since  the  fractional  rcsulua!  weight  decreases 
with  inci  casing  temperature. 

A  plot  of  {lu(  —  dn,/ilr)  —  In  f,f »■)}  versus  {1/T}  only  yields  a  straight  line  if 
f,(n*)  has  been  cliosen  correctly.  In  the  rase  or /,(»■)  —  it"  the  reaction  order  n  is 
varied  un'it  a  straight  line  is  found.  The  actis  a;ion  energy  A' can  be  obtained  from  the 
slope  of  this  plot,  the  frequency  factor  A  from  the  intercept  on  the  ordinate  and  the 
reaction  order  n  is  determined  by  the  straight  line  criterion. 

This  analj  sis  can  also  be  applied  to  the  DSC  data  since  the  basic  rqn.  (12)  only 
differs  from  (10)  by  constants. 


=  ±  All{n-0  -  w.)  A  exp  -~/T  j  f(«) 


The  -i  (  — )  in  the  equation  corresponds  to  an  exothermic  (endothermic) 
reaction  where,  Q  is  the  heal  given  off  by  the  sample,  All  is  the  enthalpy  change  of  the 

entire  reaction,  w0  is  the  initial  weight  and  w,  is  the  rest  weight,  and  the  other  para¬ 
meters  arc  as  defined  previously. 

Taking  logarithms  gives 

In  J  ~  •*'  f('v)  ~  hi  All(\vn  -  »v,)  A  -  03) 

This  time  d Qji't  at  temperature  T  is  the  measured  quantity  and  w  needs  to  be 
calculated.  This  is  easily  done  since,  if  the  reaction  results  in  a  weight  loss 


I  f  dC  ^ 

w'mAiij  dr  dl 


The  computer  numerically  integrates  dO/dr  to  yield  the  weight  at  time  I.  A  plot 
of  {In  AQ’At  -  In  f(»)}  versus  I  IT  produces  H,  A  and  n  (if  f(». )  =  O  in  the  same  way 
as  for  the  weight- loss  curves  above.  The  slope  is  again  -EjR  and  the  intercept  is 
In  d//(u  0  —  it  , 1 .4. 

In  practice  tlic  data  from  the  TG  or  DSC  is  recorded  on  8  track  ASCI  paper 
tape.  This  information  is  then  read  into  disc  files  of  live  computer.  These  files  can  be 
used  by  the  programmes  which  calculate  the  activation  energy,  etc.  After  these 
programmes  have  run,  the  plotter  output  is  viewed  on  a  television  screen.  Only  if  the 
plot  is  satisfactory  is  a  copy  obtained  from  the  plotter  for  the  final  analysis,  flic 
dimensions  correspond  to  the  original  output  of  the  chart  recorder  so  that  a  direct 
comparison  is  possible.  Different  pen  colours  arc  med  to  diflerentiate  between  the 
raw  data  (black)  and  the  theoretical  fit  to  the  dat;  (green).  These  colours  are  re¬ 
presented  in  Fig.  1  as  a  broken  line  (black)  and  n  full  line  (green). 
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fit?.  1.  Wii!;|it  toss  curve  of  C.:<  C»OjlljO  analysis.  The  rktshed  line  is  the  oiigiivd  data  and  ihefull 
line  ihc  ■  eptot  us/:;;1,  ;hc  e.  >:  p  pnteJ  (/..  A,  «)  triplet  fc*i  the  thtev  suae'.ht  line  vet  ions  \%  liiclr  arc 
svtKtratcd  hy  venical  dadw.t  tines. 


7  lie  differential  programme* 

There  urc  ihtce  different  ways  in  which  differentials  to  the  v.eisd>f-loss  curve 
were  obtained.  The  first  uses  dilfennccs  and  produces  a  noisy  differential  curve,  which 
is,  however,  neatest  to  the  original  data.  Tlte  other  two  use  a  polynomial  tit,  and  a 
spline  fit,  respectively.  The  polynomial  fa,  usually  of  25th  cider,  is  generally  pood  but 
produces  artificially  large  iimbth limits  over  straight  regions  of  the  curve.  The  best 
results  were  obtained  by  a  cubic  spline  fit.  This  piog  .•  amine  fits  a  third  order  curve  to 
sections  of  the  data,  and  joins  them  in  su-.it  way  that  the  function  and  (lie  first  two 
derivatives  are  continuous.  Tie-  programme  uses  a  standard  subroutine  from  the 
Harwell  I  orlran  library,  w.m.  h  ,  ho<<sc.\  to.  own  Vines  (points  where  two  splines  meet) 
according  to  the  b-.h;tvt<*»r  of  "m  oiviiiv;'  curve,  t  ids  means  that  over  a  straight  part 
of  the  curve  a  long  section  car:  or  titled  bv  a  single  cubic  spline,  whereas  at  places  of 
greater  variation,  smaller  sections  an  chosen  to  be  approximated.  The  differential 
of  the  curve  is  then  obtained  by  ditrcrc»:ii:«»ng  .he  third  order  spline.  The  fit  of  tlte 
spline  and  the  differentia!  is  wiy  accuiate  The  amount  of  data  smoothing  achieved 
by  this  fit  suits  the  purpose  of  the  follow  tag  programme  (next  section)  best. 

Ihc  programme  to  (bid  f  .  A  ami  n 

It  is  assumed  in  this  analysis  that  tlte  fraction  that  leads  to  the  weight-loss  curve 
is  at  least  piece-wise  of  the  nth  order,  i.c.  follows  eon.  (10)  with  f,(;r)  —  it"‘.  Ibis 
programme  plots  [In  (—  dic'd.')  -  n  In  u-}  versus  I  IT  for  different  values  of  n.  The 
values  of  dtr/Jr  and  of  u  arc  obtained  from  rt  tile  that  was  created  by  the  previous 
prof/aimne  and  contains  the  smoothed  weight -loss  data  and  the  differential  in  the 
form  of  the  spline  parameters.  The  value  of  ii  is  taken  as  zero  for  the  first  plot  and 
then  inct  eased  by  0  2  test  times  up  to  rt  -- 2 .  These  eleven  lines  have  different  curvatures. 
The  stniighicst  of  them  can  castly  be  selected  by  eye  and  this  it  value  is  then  taken  as 
the  reaction  oreci  lie'  -lope  yields  /  and  the  intercept  A.  ihc  resnitin p  triplet  of 
numbers  n)  ate  ne.essary  and  sufficient  to  characterise  an  mh-ordcr  reaction. 

The  type  of  reaction.  i  e  f(*e)  --  u*  can  be  changed  easily  so  that  the  programme  is 
applicable  to  reactions  of  any  type.  i.e.  f(v)  urn  have  any  form  (c.c.  random  chain 
scission). 

This  and  die  h'ili.mke1  iM.-giamm,-  rui  o!-t  niml  funv.  t he  .w, there 
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llw  programme  to  pm  Jure  a  theoretical  fit  horn  the  ( /.,  *f,  «;  fr/p/cf 

This  programme  allows  the  production  of  a  weight-loss  curse  with  up  to  four 
different  sections  each  of  w  hich  has  its  own  E,  A  and  n  factor.  It  was  implemented  on 
a  Hewlilt  Packard  l>;sh  top  calculator  and  plotter  of  the  tv  pc  UI’lWO?.  Tire  iheotciical 
plots  produced  in  such  a  way  fit  the  original  weight-loss  cum  very  well.  They  also 
provide  a  useful  cheek  that  all  the  calculations  necessary  to  produce  this  triplet  were 
correct,  and  produce  a  good  fit. 
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t  ig.  2.  Activation  energy  i»I©t  for  vat  ions  values  of  reaction  outer  n.  The  n  value  which  sives  i|k- 
bcsl  sliaixiil  line  M'Siucnts  is  the  otic  iipprop:iate  to  the  unction  (here  0  6V  There  me  tlnee  sups  to 
the  reaction  hut  live  fast  two  correspond  to  only  va.  2;0  of  the  ueight-ioss. 


u  st  of  mi.tiioo 

The  dehydration  reaction  of  calcium  oxalate  monohydrate 
Ca  C/)„H,0  -.  Ca  C,04  •:  H,0 

is  well  understood  and  values  for  the  reaction  parameters  can  l»e  obtained  from  the 
literature*'  It  therefore  provides  a  good  test  ease  for  the  new  analysis  method. 
Figure  1  shows  the  weight-loss  curve  of  a  sample  of  Ca  C;04H:C).  1  he  dashed  curve 
represents  the  original  data  and  the  full  curve  the  theoretical  fit  based  on  the  computed 
(E,  A,  n)  triplet.  Figure  2  shows  the  activation  enerpy  plot  The  number  at  the  cm* 
each  line  gives  the  reaction  order  that  is  assumed  for  the  ca  leu  It:  lion  of  the  emu 
sponding  line.  The  straight  lines  which  approximate  the  cur  ves  best  arc  used  to  obtain 
(£,  A,  n)  triplets  for  the  three  sections.  Since  these  (E.  .1.  n)  triplets  completely 
specify  the  reaction  the  original  wcight-ljss curve  can  be  modelled  (full  curve  fig.  1 V 
The  results  obtained  with  this  method,  and  literature  values  arc  summarized  in  Table 
I.  All  our  experiments were  performed  at  a  heating  rate  of  10  den  min" '.  The  results 
obtained  in  this  way  have  an  error  of  -  5  kJ  mol ' 1  w  hich  is  due  to  the  variance  of  the 
samples  themselves.  The  method  itself  has  much  greater  accuracy.  It  was  found  that 
the  best  straight  line  could  be  fitted  to  a  curve  corresponding  to  a  reaction  order  of  0  6 

TABLE  1 


E  l/i  4  Ref. 

(kJ  mole-'} 


*7  19.8  Ibis  work 

84  I7.S  This  work 

88  21.0  * 

97  9 
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which  indicates  a  surface  reaction  (theoretically  2/3).  It  can  also  he  seen  on  Fig.  2 
(hat  there  arc  three  distinct  stages  of  the  reaction.  The  first  two,  however,  correspond 
to  only  2%  of  the  weight  and  are  attributed  to  volatile  substances  on  the  surface. 
This  shows  the  sensitivity  of  the  method  used  especially  considering  that  the  vertical 
axis  has  a  logarithmic  scale.  The  real  dehydration  reaction  is  represented  by  the  Iasi 
straight  line  segment. 

The  method  outlined  in  this  paper  has  recently  been  successfully  applied10  to 
decomposition  studies  of  P.E.T.N.,  P.F..T.N.  with  fillers,  high  density  polyethylene 
(HOPE)  and  polylctrafluorocthylcnc  (PTFL). 
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SECTION  4 


ANALYSIS  OF  ISOTHERMAL  KINETIC  DATA 
FROM  SOLID-STATE  REACTIONS 

Tong  B.  Tanc;  and  M.  M.  Ciimjdhki 

PhysUi  ami  C'/iemistry  of  Solids,  Cucrmlish  Laboratory,  Mailing!*  y  Road,  Can, hi  ijge  UK 
(Received  December  12.  IM78) 


In  most  solid  state  .  tactions  the  reaction  velocity  can  be  described  as  a  product  of 
two  functions  K(i  )  an  J  /<1  -  1)  where  7  is  the  temperature  and  .  the  degree  of  con¬ 
version  of  the  solid  reactant,  t  he  physical  interpretation  of  these  functions  is  discussed, 
and  a  systematic  met  nod  is  described  by  winch  A I  A  of  a  reaction  i»  identified  from 
its  kinetic  data.  A(7")  and  the  icaction  mechanism  arc  then  determined.  *1  hi  '  . -d 

has  been  successfully  applied  to  analyse  the  kinetics  of  the  thci mat  dcconipo.it.  n  of 
silver  azide. 

In  a  solid-state  rcactmn,  the  reaction  velocity  is  given  by  -  d/d;  (I  —  a)  *  a, 
where  a  =  a(f)  is  the  fraction  of  the  solid  roaclant  winch  has  reacted  by  time  /. 
Its  kinetics  can  be  solved  by  determining  i  as  a  function  of  temperature  and  the 
global  amount  of  reactant  left.  This  phenomenological  knowledge  is  a  necessary, 
though  not  sufficient,  condition  f<»i  elucidating  the  reaction  mechanism.  It  is 
necessary  to  formulate  the  reaction  velocity  in  terms  of  the  global  variable  a, 
because  iherc  is  a  continuous  collapse  of  structure  in  ilie  reactant  Furthermore, 
the  local  concentration  of  '.he  reactant  varies  throughout  the  reaction  volume  and 
cannot  be  used  as  a  state  variable.  In  fact,  unlike  ihc  case  of  a  homogeneous  re¬ 
action  in  the  liquid  or  gaseous  phase,  there  is  no  real  ‘reaction  order’  with  respect 
to  any  reactant  in  a  reaction  involving  condensed  matter  whose  mechanism  is 
usually  of  the  heterogeneous  type. 

If  the  reaction  proceeds  isothcrmally,  it  is  observed  empirically  iltat  (a,  t) 
curves  corresponding  to  different  temperatures  T  arc  isomorphic  to  one  another, 
at  least  within  a  range  of  T,  i.c..  by  a  linear  scale  change  in  different  curves  can  be 
superimposed  [1).  It  follows  that  i  is  a  separable  function: 

~  ~  2).  (I) 

Here /(I  -  a)  may  change  in  different  ranges  of  T  or  a.  For  every  /(I  -  a),  there 
corresponds  a  single  A'(’T).  It  should  be  noted  that  experimental  data  may  be  ade¬ 
quately  analysed  by  U)  only  if  it  has  been  ensured  that  the  tempet  aturc  distribution 
in  the  sample  is  sufficiently  uniform  anti  constant.  Furthermore,  the  theoretical 
significance  of  K(T)  and  /(I  —  7) determined  lroi  1  the  data  should  always  be  exam¬ 
ined  withj  tegard  toi  the  class  of  mechanisms  .hey  indicate.  A  question  of  con- 
istency  arises  in  this  respect.  In  the  literature,  coherent  and  integrated  accounts 
of  the  physical  (as  contrasted  with  the  formal)  meanings  of  A'( T)  and  /(I  -  2) 
are  not  easily  found.  A  discussion  on  their  ir.erprctations  therefore  forms  the 
first  part  of  this  paper.  In  the  second  part  we  propose  an  efficient  method  of  de¬ 
termining.  with  confidence,  both  A I  -  2)  and  K  T)  from  y  and  2(1)  data.  The  inter¬ 
pretation  of  JO  —  2)  is  the  essential  link  in  this  method  between  the  experimental 
data  and  the  functions. 


Physical  interprets  ions 


The  function  f(  I  —y  I 


Solid-state  reactions  are  complex  process...  which  proceed  in  scvctal  stages. 
These  can  be  Ihc  dvlocalnation  or  u.m  .lei  o'  ,  lections  in  chemical  bonds  (it.  the 
ease  of  non-metals),  the  dilftision  of  atoms,  lie  radicals,  or  ions,  the  desorption 
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of  product  molecules  when  they  are  in  the  gaseous  phase,  the  heat  transfer  to  the 
reaction  zone  in  the  case  of  endothermic  reactions,  and  the  formation  of  a  new  solid 
structure  (crystalline  or  amorphous)  if  one  of  the  products  is  in  the  solid  phase. 
The  last  step  may  often  be  further  differentiated  into  nucleation.  growth  of  nuclei 
(at  velocities  which  depend  on  sizes  of  nuclei)  [I]  and  sometimes  the  collapse  of 
the  lattice  from  a  transitory  one  to  the  equilibrium  structure  [2].  Irrespective 
of  the  details  of  reaction  mechanisms,  however,  under  a  given  set  of  circumstances 
(T.  x.  sample  history  etc.)  one  of  the  stages  will  be  the  slowest.  It  then  acts  as  the 
rate  limiting  step  of  the  reaction,  and  it  will  determine  the  kinetics  i.e..  the  rate 
law  in  ( I ). 

Further,  a  solid-state  reaction  has,  in  contrast  to  a  homogeneous  reaction  whose 
progress  is  independent  of  spatial  coordinates,  an  additional  controlling  factor, 
namely  topochemistry.  This  refers  to  the  geometrical  shape  of  the  solid  reactant 
and.  in  different  cases,  to  its  free  surface  area,  its  defect  structure,  the  thickness 
of  the  product  layer  if  solid,  or  to  the  product- reactant  boundary,  etc. 

The  function  A 1  —  *)  reflects  the  nature  of  the  rate-limiting  step  and  the  topo¬ 
chemistry  of  the  reaction.  It  may  accordingly  depend  on  certain  sample  conditions, 
such  as  whether  the  sample  has  been  pre-irradiated  or  bleached,  and  whether 
the  sample  is  in  the  form  of  a  powder,  or  a  large  single-crystal  of  a  different  shape 
from  the  crystallites.  It  will  vary  in  several  ranges  of  T  if  in  each  of  them  a  different 
elementary  step  becomes  rate-limiting,  as  occurs  in  the  decomposition  of  potassium 
azide  [3].  At  a  given  temperature,  it  may  also  change  in  different  ranges  of  a. 
due  to  the  switching  of  the  rate-limiting  step  or  topochemical  changes.  This  hap¬ 
pens.  for  instance,  in  the  oxidation  of  zirconium  [4]  and  in  most  decomposition 
processes  [1],  An  extreme  case  is  the  decomposition  of  ammonium  perchlorate 
which,  in  the  two  temperature  regimes  below  and  above  620K,  has  entirely  dif¬ 
ferent  reaction  mechanisms  and  in  fact  yields  different  reaction  products  [5], 
In  such  cases  there  may  be  competing  paths  for  the  chemical  reaction  or  it  may 
in  fact  be  followed  by  another  chemical  reaction  whose  ‘onset  temperature’  is 
higher.  In  all  possibilities,  however, /Tl  -  a)  should  be  the  same  for  a  given  a 

independent  of  /'  (within  a  range)  if  it  is  to  have  more  than  only  an  empirical 
significance. 

In  Table  1  we  have  collected  together  the  more  common  forms  of  ft  I  -  a) 
which  have  been  used  in  the  literature,  and  the  corresponding  integrated  forms 

i  dr  iifil  -  a)  «  j  da//(i  -  a)  s  F(a).  Note  that  F( a)  =  K  (r  -  i„)  if  the  range 
o  6 

of  a  for  which  it  becomes  applicable  starts  at  a  *  x(rn),  Also,  for  simplicity  here¬ 
after  we  write  K  for  K(T). 

In  many  reactions,  such  as  most  deccn-pcsirirrs  and  dehydrations,  the  rate- 
limiting  step  takes  place  at  the  interface  tctwccn  different  phases  as  in  sublimation. 
The  speed  at  which  the  interface  imvi  UiC  *  s actant  is  (at  a  given  temperature) 
then  either  a  constant,  or  a  unique  function  of  the  interfacial  area.  This  area 
therefore,  from  the  kinetic  point  of  view,  plays  the  same  role  as  that  of  con¬ 
centration  in  homogeneous  reactions.  If  the  speed  is  constant,  then  the  theoretical 
significance  of/U  -  a)  is  clear:  it  gives  the  area  expressed  as  a  fraction  of  the 
original  area  at  *  =  0.  This  is  the  case  of  a  reaction  controlled  by  the  movement 
of  a  coherent  phase-boundary  and  listed  as  F.  G  and  H  in  Table  1.  In  this  situation, 
the  explicit  form  of /(I  -  a)  depends  on  the  geometry  of  the  reacting  system, 
though  generally  it  is  a  decreasing  function  of  a  or  at  most  constant. 

If  the  reaction  consists  of  the  formation  of  compact  nuclei  of  a  solid  product 
at  localized  places  in  the  reactant  followed  by  their  relatively  rapid  growth,  then, 
to  express  the  total  interfacial  area./fl  -  a)  is  derived  from  the  laws  of  nucleation 
and  growth.  This  is  the  situation  when  the  reaction  is  autocatalytic  (13]:  reactant 
molecules  at  a  reactant-product  interface  react  in  preference  to  those  at  a  reactant- 
‘  vacuum'  surface.  The  preference  is  due  to  the  existence  of  microstrains  in  the 
reactant  at  the  interface,  or  due  to  the  electrochemical  potential  of  the  product 
phase  when  the  rate-limiting  step  is  a  redox  process.  The  various  possible  forms  of 
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/(I  —  x)  for  an  autocatahtic  reaction  are  listed  in  Table  I  .A  to  £.  Note  that  they 
have  the  general  form  /(I  -  at)  =  »p(l  -  x)q.  In  A.  B  and  £,  q  is  zero  and  i 
increases  monotonically  with  a.  Such  a  situation  is  most  unlikely  to  last  up  to 
z=l.  These  types  of  /(I  -  a)  therefore  may  apply  only  to  the  acceleratory  part, 
if  it  is  present  in  the  a-time  curve  and  which  will  usually  be  followed  by  a  decay 
part.  In  C  and  D.  q  is  non-zero  and  these  types  give  sigmoid-shaped  curves.  The 
inflexion  point  occurs  at  a'  =  p(p  +  q).  as  can  be  seen  at  once  from  the  condition 
i  =  0.  It  is  thus  kinetically  feasible  for  them  to  fit  the  complete  experimental  curve. 
There  are.  however,  physical  grounds  to  consider  that  even  they  should  be  used 
to  analyze  only  the  acceleratory  period  [I]. 

In  other  reactions  the  rate-limiting  step  is  not  confined  to.  or  does  not  only 
occur  at.  the  reactant  surface.  For  instance,  in  the  unimolccular-decay  type  of 
reaction,  all  molecules  whether  on  the  surface  or  in  the  bulk  have  an  equal  proba¬ 
bility  per  unit  time  of  reacting.  This  is  the  case  when  the  change  from  the  re¬ 
actant  to  the  solid  product  phase  involves  little  re-arrangement  of  the  reactant 
atoms.  The  reaction  has  homogeneous  mechanism  and  thus  a  true  reaction 
order  of  one.  i.e.  f{  I  -  x)  =  1  —  x.  Many  decomposition  reactions  tend  to  this 

Table  1 

The  common  types  of  solid-state  reaction 

Reaction 


Autocatalytic  Power-law  nucleation  and  growth  at  constant  speed:  A 

Linear  branching  chain  of  nuclei,  no  overlap  during  growth:  B 
Branching-chain  nucleation.  interference  during  growth:  C 
Random  nucleation,  growth  accompanied  by  ingestion  of  nuclei:  D 
Instantaneous  nucleation.  size-dependent  growth:  E 

Phase- boundary  controlled  decay  1  -dimensional:  F 

2- dimensional:  C 

3- dimensional  H 

Unimolecular  decay :  / 

1 -dimensional:  J 

Diffusion-controlled  2-dimensional:  K 

3-dimensional:  L 


At  -  «J  -  i:K 

r is>  -  Kt 

Reference 
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0.57 
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i)  Log  J  — In  (I  -  »)}  Constant  +  m  log  f,  0.15  <  x  <  0.5 

ii)  Plot  of  L.H.S.  in  if  against  log  t  distinctively  concave  upwards:  C„  and  x„  arc  constants 
while  C(T)  is  function  of  temperature 

iii)  b  =  0.774.  0.700.  0.664.  0.642 .  0.556  for  n  =  2.  3.  4.  5 _ x 

is)  Alternative  derivations;  we  have  obtained  the  approximate  forms  of  fit  -  x)  by  ex¬ 
panding  into  series  (I  —  r)-5'*  and  (I  -  z)_M  to  second  order  in  »:  resulting  error  at  x*  6 
(<10*;  for  j  <  0.*> 
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limit  at  high  a  values.  Another  example  is  when  the  rate-limiting  step  is  the 
migration  of  product  ions  along  the  dislocation  network  to  form  additional 
growth  nuclei  at  dislocation  nodes  [17].  This  leads  to /[l  -  a)  =  %,  the  same  form 
as  for  branching  nuclei  [7].  The  slow  process  of ‘ageing’  in  some  explosives  when 
they  are  stored  at  room  temperature  may  be  by  such  a  nucleus-chain  mechanism. 
A  third  category  consists  of  reactions  controlled  by  the  diffusion  of  reactants 
across  a  product  layer  which  is  solid.  The  diffusion  may  proceed  uniformly  through 
the  bulk  of  the  layer  and  is  thus  structure-insensitive,  or  preferentially  along  its 
gross  lattice  imperfections  arising  from  product-reactant  mismatch.  In  the  case 
of  uniform  diffusion,  the  speed  at  which  the  product-reactant  interface  moves  is 
a  function  only  of  the  product  thickness  (and  temperature),  and  the  appropriate 
forms  of/(l  —  a)  are  included  in  Table  1  as  J,  K  and  L.  The  oxidation  of  metals 
often  follows  diffusion-controlled  kinetics;  in  sheet  form  these  tarnish  according 
to  the  parabolic  law  or  x.t.  Exceptions  aret  hose  metals  in  Groups  la  and  1  la  of 
the  Periodic  Table.  Excluding  beryllium,  they  all  form  oxide  layers  which  are 
porous,  so  that  the  atoms  of  the  metal  do  not  have  to  diffuse  through  a  coherent 
layer  before  coming  into  contact  with  oxygen. 


The  function  K(T) 

It  is  almost  always  the  case  that  the  temperature-dependent  part  of  (1)  can  be 
represented  successfully  by :  - 

K{T)  =  K,  exp  (  —  EjkT)  (2) 

in  which  k  is  Boltzmann’s  constant,  and  the  macro-kinetic  constants  E  and  Kr_ 
do  not  depend  on  T  (within  the  range),  though  usually  they  take  on  different  values 
when  /( 1  -  a)  changes. 

If  it  is  established  that  the  reaction  is  rate-limited  by  a  diffusion  or  migration 
process,  the  interpretation  of  K(T)  is  complicated,  but  obviously  it  is  proportional 
to  the  corresponding  transporfcoefficient,  which  in  general  is  an  exponential  func¬ 
tion  of  T.  An  over-simplified  theory  for  the  situation  of  uniform  one-dimensional 
diffusion  gives  K(T)  =  (SI  Vn)2D(T),  where  S  is  the  interfacial  area.  V„  the  initial 
volume  of  the  reactant,  and  D(T)  the  diffusion  coefficient  (cf.  [19]). 

For  a  single-solid-reactant  reaction  which  is  controlled  by  a  surface  process, 
on  the  other  hand,  the  simple  theory  of  Shannon  [18]  is  often  successful.  This 
theory  is  a  generalization  of  the  Polanyi  —  Wigner  equation.  Assuming  the  existence 
of  some  activated  complex,  which  as  a  transition  state  can  be  treated  in  thermo¬ 
dynamic  equilibrium  with  the  reactant,  he  related  the  pre-exponential  factor  K 
to  the  rotational  and  other  internal  degrees  of  freedom  of  a  reactant  molecule 
in  addition  to  the  vibrational  ones.  Following  Shannon  we  can  set;  - 


=  (kT/h)  exp  (ASx/k)  dSJV0 

(3) 

and 

exp(  —  E/kT)  =  exp  (—  AH:/kT). 

(4) 

Here  the  mean-frequency  factor  kT/h  containing  Planck's  constant  is  usually  in 
the  region  of  lO^s^fsee  below).  dS*  and  AHX  are  respectively  the  entropy  and  the 
enthalpy  of  formation  of  the  transition  complex,  d  is  the  thickness  of  one  mono- 
layer  and  Vn  the  initial  volume  of  the  reactant,  and  S/(l  -  a)  gives  the  free  sur¬ 
face  or  the  product-reactant  interface  area  when  the  degree  of  conversion  is  a. 
(Strictly  speaking,  it  has  been  assumed  that  the  reaction  proceeds  isobarically). 

Note  that  in  this  interpretation  the  empirical  quantity  Kr  contains  the  surface- 
-to-volume  ratio  and  so  depends  on  the  sample  geometry.  Also,  it  is  apparently 
proportional  to  T.  (In  gas  reactions,  the  collision  theory  gives  Kr  xT1.)  In  our 
opinion,  however,  if  the  vibrational  modes  are  being  considered  then  only  at  low 
temperatures  will  the  peak  distribution  of  phonon  frequencies  lie  at  kT/h.  For 
most  substances  (with  the  exceptions  of  Be.  Cr  and  diamond)  the  Debye  tempera¬ 
ture  fl0  is  less  than  500K.  so  that  the  frequency  factor  should  stay  as  kOn  li  <  lO'^s* 1 
for  all  likely  experimental  temperatures. 
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The  factor  exp(dS;  k)  may  alternatively  be  written  in  terms  of  partition  func¬ 
tions  as  Q*!Q .  which  can  be  determined  from  spectroscopic  data  [20].  In  most 
cases  dS5  cannot  be  larger  than  the  reactant  entropy  of  melting,  and  exp(ASt/k) 
comes  out  normally  between  unity  and  104.  Occasionally  exp(dS*/A)  is  found  to 
be  less  than  unity,  as  is  the  steric  factor  in  gas  reactions.  Such  a  negative  value  of 
JS*  means  that  the  activated  complex  is  more  ordered  than  the  reactant  (e.g. 
[2 1  ]).  Experiments  on  some  decomposition  reactions  have  given  K  which  are  ab¬ 
normally  high  in  comparison  to  the  theoretical  values  of  (3).  Hypotheses  put  for¬ 
ward  to  explain  such  discrepancies  include  co-operative  activation  [22],  proton- 
delocalization  [23].  and  a  mobile  layer  of  molecules  on  the  reactant  surface  [18]. 


Kinetic  analysis 

Current  practice 

In  determining  the  kinetics  one  wishes  to  find  K , .  E.  and  /( 1  —  at)  or  equiv¬ 
alently  F( a)  so  that  the  reaction  velocity  can  be  predicted  at  any  given  T  and  a. 
This  is  commonly  done  by  analyzing  a  set  of  ae<r )  or  equivalently  a(/)  values  ob¬ 
tained  by  monitoring  a  number  of  samples  reacting  isothermally  at  a  number  of 
temperatures.  The  consistency  of  the  K *  and  E  values  with /( 1  —  a)  should  as  far 
as  possible  be  assessed,  and  correlated  with,  for  instance,  microscopy. 

A  quick  method  of  calculating  E  was  used  by  Haynes  and  Young  [24].  Consider 
a  set  of  (a.  t)  curves  which  have  been  found  to  be  isomorphic.  For  any  two  curves 
(a,,  ty)  and  (22,  /*)  corresponding  to  temperatures  7",  and  T.t  respectively,  one  can 
write 

F{ a,)  =  tyK,  exp  (-E/AT,) 

F(*4)  =  t,Kx  exp  (-EAT,). 

By  choosing  points  corresponding  to  the  same  a  on  the  two  curves  so  that  F{ a,)  = 
=  E(a2).  £  can  be  evaluated  by  plotting  In  /  vs.  1  T.  On  the  other  hand,  to  deter¬ 
mine  F(ot)  often  a  trial-and -error  method  is  resorted  to.  Conflicting  forms  of  the 
function  have  sometimes  been  asserted  by  several  authors  for  the  same  material, 
like  NH4C10<  (see  [25])  and  KMnO,  (see  [26]). 

A  conventional  way  of  superimposing  isothermal  curves  is  to  convert  them  into 
‘reduced-time  plots'  by  individually  scaling  their  r-axis  with  the  factor  l/r(.  where 
tj  is  the  time  when  *  =  0.5  on  the  ;-th  curve.  In  this  way  K{T,)  is  absorbed  into 
each  scale  factor  and  all  (s,  r)  =  (0.5.  1)  points  coalesce,  while  other  *(/)  points 
may  be  plotted  out  to  see  if  the  curves  are  indeed  isomorphic.  Sharp  et  al.  [27] 
tabulate  the  theoretical  values  of  a  against  t'x  for  some  of  the  E(a)  shown  in  Table  1. 
They  propose  that  by  comparing  experimental  data  with  such  master  values  the 
correct  F(oc)  can  be  identified. 

The  above  method  may,  however,  result  in  ambiguity  due  to  a  number  of  as¬ 
pects.  Experimental  data  contain  random  errors,  but  no  simple  statistical  analysis 
can  be  applied  to  the  identification  criterion  it  employs  because  no  straight-line 
graphs  are  involved.  Additionally,  a  general  problem  for  all  isothermal  experi¬ 
ments  is  the  zero-time  uncertainty.  The  finite  time  taken  by  the  sample  to  reach 
the  designated  temperature  may  be  negligible  relative  to  r.  yet  may  affect  the  com¬ 
parison  with  the  tabulated  values  [28].  Moreover.  FM  may  change  in  different 
regimes  of  the  a (r)  curves,  as  mentioned  earlier. 


A  new  method 

Here  we  suggest  a  step-by-step  approach  to  determine  F{x).  It  was  noted  by 
Hancock  and  Sharp  [28]  that  for  many  forms  of  £(*».  the  plot  of  log[-ln(I  -  *)] 
vs.  log  f  is  almost  linear  if  *  is  restricted  to  between  0.15  and  0.5.  Using  a  com¬ 
puter  program  to  generate  artificial  values  and  their  log-ln  plots,  we  have  found 
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this  true  for  all  the  theoretical  forms  listed  in  Table  1.  with  the  exception  of  B, 
C  and  E.  The  slope  in  each  case  is  listed  there  under  the  Column  'in' . 

Obviously  a  log-ln  plot  is  not  very  sensitive.  If  we  were  to  rely  solely  on  it  to 
discriminate  between  the  functional  forms  of  £(x).  the  experimental  data  would 
have  to  be  of  the  highest  quality.  A  slightly  more  sensitive  way  is  to  plot  [-  ln(  1  — 
-  it]' is.  i.  but  then  m  can  only  be  obtained  iteratively.  Fortunately,  m  does 
differ  significantly  between  different  groups  of  A(x).  and  the  final  discrimination 
is  easily  achieved  by  a  further  graphical  step.  There  are  three  possible  situations 
for  this  second  step:  - 

in  ’  or  log-ln  concave  upwards 

It  will  be  seen  from  Table  I  that  this  situation  suggests  an  autocatalytic  reaction 
tor  which  x  -  Axr(l  -  xt*1  for  certain  />  and  q.  From  the  experimental  data  of  x 
and  /  one  can  then  do  a  least-squares  fit  on  the  graph  of  f  log  -/  -f  log  x  against 
.1  logi  I  -  x)  .1  log  x.  and  find  />  from  the  r-intercepl  and  q  from  the  slope.  Here 
I  log  x  =  log  x(/,)  -  log  x(M.  etc. 

For  an  x  expression  of  this  form,  one  has  />(  I  -  y.  )  =  qy  .  where  x'  is  the  value 
at  maximum  x.  Using  this  relation  to  reduce  the  number  of  unknown  parameters 
to  one.  one  can  use  a  simpler  graph  to  determine  p  and  q  [26].  However,  the  cal¬ 
culated  values  of  p  and  q  are  then  subject  to  the  accuracy  of  x'  and.  more  funda¬ 
mentally.  the  possibility  that  p  and  q  may  change  from  one  range  of  x  to  another 
is  not  allowed  for.  As  mentioned  above,  those  types  of  f{\  -  x)  in  which  q  —  0 
represent  the  accelerators’  period  which,  in  general  should  be  followed  by  a  decay 
period  governed  by  a  different  form  of  /i  I  -  x). 


The  reaction  is  either  phase-boundary  controlled  or  ummolecular.  and  x  = 
=  A'(l  -  x)\  as  seen  in  Table  I.  One  then  draws  the  graph  of  log  x  against 
log(l  -  x)  to  find  r.  the  apparent  reaction  order. 

3.  m  -  0.5 

The  reaction  is  diffusion  controlled  (see  Table  1).  One  has  to  test  separately 
whether  xU)  is  parabolic  (the  diffusion  is  in  one  dimension),  or  —  x  ln(l  -  x)  = 
=  K(  I  —  xf  with  s  =  0  (two  dimensions)  or  s  =  1  3  (three  dimensions). 

The  correlation  coefficient  in  the  least-squares  fit  serves,  by  measuring  the  li¬ 
nearity  of  the  x  graph,  as  a  quantitative  indication  of  the  confidence  to  be  attached 
to  the  identified  form  of /(I  -  xl.  It  may  be  that  p  and  q.  r.  or  s  change  once  or 
twice  as  the  reaction  proceeds  from  beginning  to  completion,  but  the  i  graphs 
will  show  it  by  displaying  several  linear  segments.  If  however,  a  part  of  the  graph 
say  from  (x,.  t,)  to  (x.,.  t.,)  is  non-linear,  then  F(x)  has  changed  to  a  form  in  another 
group.  The  first  step  should  then  be  repeated  for  that  part:  log  [—  Inf  I  -  [x  +  x,  ] ' 
[I  -  x,])]  is  plotted  rs.  !og(r  -  /,)  for  x  between  x,  -+-  0. 1 5(x_.  -  x,)  and  x,  +  0.5 
(x4  -  x,).  followed  by  one  of  the  above  three  alternative  procedures.  On  the  other 
hand,  if  a  good  fit  is  found  with  values  of  p  and  q.  r.  or  s  that  are  not  in  Table  1. 
the  experimenter  should  assess  whether  theoretical  justification  can  be  provided. 
In  this  way  new  rate  laws  may  be  identified. 

Confirmation  is  carried  out  by  plotting  the  selected  functional  form  or  forms 
on  top  of  the  experimental  curve.  A  slight  misfit  in  the  very  early  part  (-  1 
minute,  depending  on  the  sample  size  and  the  environment)  may  be  attributed 
to  thermal  lag-time  and  ignored.  The  x  graphs  give  A',  and  from  a  set  of  A  values 
at  different  temperatures  A  ,  and  £  can  be  determined.  Note  that  the  determination 
of  these  macro-kinetic  constants  depends  on  the  form  of  F(l  -  x)  chosen,  as  it 
should  be. 
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(t  may  be  added  that  we  see  it  an  immediate  possibility  to  have  full  automation 
in  the  acquisition  and  processing  of  data  in  thermal  analysis  experiments.  The 
hardware  can  be  under  the  control  of  microprocessors  or  dedicated  minicomputers, 
and  their  output  would  go  into  a  computer  or  the  same  minicomputer.  A  com¬ 
puter  program  can  then  reduce  the  data  to  x(t)  or  Hr)  curves,  and  further  analyse 
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Fig.  I.  Decomposition  of  a  single  crystal  of  AgN,  at  551  K.;  1.  Thermogravimetric  data. 
2.  jr-timc  curve.  3.  1  -  (1  —  *)' :  is.  time. 
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Fig.  2.  Log-In  plot 


the  curves  to  identify/!;  1  -  a)  and  so  calculate  Kr  and  £.  according  to  the  method 
proposed  here.  Nevertheless,  the  physical  interpretation  of  these  results  by  the 
experimenter  remains  the  crucial  step. 

The  method  has  been  applied  to  investigate  the  kinetics  of  slow  thermal  de¬ 
composition  in  silver  azide  single  crystals  using  thermogravimetric  data.  Curve  ! 
in  Fig.  1  is  a  typical  experimental  curve  of  weight  loss  against  time  t,  and  Curve  2 
is  the  corresponding  reduced-time  plot  of  x  rs.  r  t.  In  Fig.  2  we  plot  log[-ln(l  - 
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Fig.  V  y  graph 


—  2)]  vs.  log ((/t)  for  0.15  <  2  <  0.5.  obtaining  essentially  a  linear  graph  of  slope 
~1.1.  The  decomposition  of  AgN:l  therefore  appears  to  be  phase-boundary 
controlled  or  unimolecular.  Accordingly,  we  draw  in  Fig.  3  the  graph  of  log  2  vs. 
log C l  —  2).  It  shows  that  between  1—2=1  to  0.1.  z  =  All  —  z)-  with  the 
correlation  coefficient  among  the  data  points  better  than  +0.9.  Indeed,  the  plot 
of  1  —  (1  —  z)!  vs.  t  (Curve  3.  Fig.  1)  is  a  good  straight  line,  with  a  correlation 
coefficient  of  +0.99  for  0.1  <  2  <  0.9.  The  decomposition  is  thus  of  the  contract¬ 
ing  cylinder  type.  We  have  also  obtained  K ,  and  E.  and  found  that  their  inter¬ 
pretation  in  terms  of  Eqs  (31  and  (4)  leads  to  a  plausible  physical  picture.  Further 
details  are  given  in  a  separate  paper  devoted  to  the  mechanism  of  the  thermal  de¬ 
composition  of  AgN:,  [30], 

In  a  second  paper  [29],  we  discuss  the  analysis  of  dynamic  kinetic  data  which 
are  easily  obtained  from  modern  thermoanalvtical  instruments. 


Conclusion 

We  have  described  a  systematic  method  for  determining  the  kinetics  of  solid- 
state  reactions  from  isothermal  data.  It  may  be  stressed  once  more  that  whenever 
possible  a  judgement  should  then  be  made  on  the  consistency  of  the  gnd  £ 
values  with  the  implication  of  /(l  -  2)  regarding  the  likely  mechanism  of  the 
reaction. 
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SECTION  5 


ANALYSIS  OF  DYNAMIC  KINETIC  DATA  FROM  SOLID-STATE 

REACTIONS 
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(Received  December  12,  1978;  in  revised  form  August  20,  1979) 


The  kinetics  of  heterogeneous  reactions,  involving  one  reactant  in  the  solid  phase, 
usually  follow  the  law  x  =  Km  exp(—  EjkTrfi  1  —  x),  where  x  is  the  degree  of  conversion 
of  the  solid,  and  Kw  and  E  are  the  kinetic  constants.  A  critical  examination  is  given  of 
the  various  methods  which  are  currently  used  to  analyse  dynamic  experimental  data. 
The  limitations  of  these  methods  and  their  insensitivity  to  the  form  of  /(I  —  x)  are 
pointed  out.  An  alternative  approach  free  from  these  limitations  is  suggested.  In  this, 
/(I  —  x)  is  determined  from  isothermal  experiments,  and  then  the  dynamic  data  are 
accurately  analyzed  to  obtain  the  values  of  the  kinetic  constants.  A  case  study  is  given 
to  elucidate  the  applicability  of  the  approach. 

There  are  many  reactions  of  interest  in  which  one  of  the  reactants  is  in  the  solid 
phase.  These  reactions  can  be  classified  variously  as  decomposition,  dehydration, 
calcination,  dehydroxylation,  reduction,  polymeric  inversion  and  degradation, 
oxidation  etc.,  and  they  occur  in  a  wide  range  of  substances  including  ceramics, 
explosives  and  biological  materials.  The  dynamic  method  of  studying  their  kinetics 
involves  measuring  the  reaction  rates  under  conditions  of  a  continuous  temperature 
change. 

Originally  proposed  by  Skramovsky  [1],  the  dynamic  method  is  becoming  in¬ 
creasingly  popular,  especially  with  the  development  of  differential  thermoanalytical 
techniques  like  DSC  and  DTA.  Since  the  initial  temperature  can  be  chosen  so  that 
the  reaction  rate  is  relatively  insignificant  to  begin  with,  it  does  not  suffer  from 
zero-time  inaccuracy -a  problem  which  exists  in  isothermal  experiments  where  the 
temperature  is  raised  rapidly  and  then  held  constant  at  a  particular  value.  A  further 
advantage  is  that,  provided  the  dynamic  data  have  been  unambiguously  and  cor¬ 
rectly  analyzed,  any  changes  in  the  kinetic  constants  will  not  be  over-looked  even 
within  small  temperature  intervals.  In  contrast,  the  isothermal  method  only  pro¬ 
vides  values  averaged  over  discrete  points  in  temperature.  Also,  when  the  method  of 
analysis  used  is  such  that  the  kinetic  constants  are  calculated  from  each  dynamic 
curve  then  very  few  samples  are  required ;  only  a  milligram  or  so  of  the  material  is 


needed  for  its  thermal  characterization.  If  many  runs  are  indeed  carried  out.  differ¬ 
ences  between  individual  samples  can  be  determined.  The  last  two  advantages  are 
particularly  useful  in  single  crystal  work. 

On  the  other  hand,  intrinsic  differences  should  be  carefully  distinguished  from 
the  effects  of  experimental  conditions.  In  the  first  case  literature  data  have  shown 

that  experimental  parameters  such  as  sample  mass  [53]  and  shape,  particle  size  in 
the  case  of  powder  samples,  and  ambient  atmospherefor  vacuumlcan  affect  signif¬ 
icantly  the  calculated  values  of  kinetic  constants  (see  [2]).  If  this  happens,  then 
whenever  possible  the  empirical  results  should  be  extrapolated  to  refer  to  a  stan¬ 
dard'  set  of  experimental  conditions.  Secondly,  the  healing  rate  very  often  affects 
the  shape  of  the  dynamic  curve  obtained,  but  discussion  will  be  deferred  to  the  next 
section.  Lastly,  the  very  fact  that  temperature  is  now  a  variable,  in  addition  to  time, 
complicates  the  analysis  of  data.  If  due  care  is  not  taken,  either  inaccurate  or  totally 
misleading  values  are  obtained.  In  fact,  a  survey  of  the  literature  reveals  several 
instances  of  high-quality  experimental  data  being  mis-interpreted  by  methods 
beyond  their  ranges  of  validity.  In  this  paper,  we  first  describe  the  various  methods 
and  point  out  their  limitations.  We  then  defend  the  approach  in  which  use  is  made 
of  both  isothermal  and  dynamic  experiments.  The  analysis  of  isothermal  data  yields 
/(I  -x)  (as  defined  below)  unambiguously  and.  knowing/!  1  -  a),  one  can  calculate 
individual  values  of  the  kinetic  constants  from  each  set  of  dynamic  data.  We  may 
mention  that,  historically,  isothermal  experiments  were  the  only  ones  employed 
in  the  pioneering  age  in  the  twenties  and  thirties,  when  solid-state  reactions  began 
to  be  studied  from  the  modern  point  of  view,  as  distinct  from  that  of  Langmuir. 
Nernst  and  Tammann. 

Kinetic  equation 

As  discussed  elsewhere  [3]  the  kinetics  of  a  reaction  proceeding  isothermally  can 
usually  be  described  by  the  empirical  relation : 

(“Uth.rm.1  =  ^r/d  -a)exp(-E/kD.  (IV 

Here  i  is  the  fraction  of  the  solid  reacted,  k  Boltzmann's  constant.  T  the  tempera¬ 
ture.  and  /(l  -x)  and  the  constants  K,  and  £  are  characteristic  to  the  reaction. 
The  function  /( 1  -  x)  may  change  in  different  ranges  of  x  but  is.  for  a  given  x,  inde¬ 
pendent  of  T.  at  least  within  a  range  of  T.  K-,  and  £  should  be  the  same  for  the 
same /( 1  -  x).  If  the  rate-controlling  step  of  the  reaction  occurs  on  the  reactant-free 
surface  or  on  the  reactant-product  (solid)  interface,  then  Kr  will  contain  the  sur- 
face-to-volume  ratio.  In  other  words,  the  reacting  system  should  really  be  normal¬ 
ized  per  unit  area  rather  than  per  unit  size,  and  Kr  be  given  in  units  such  as  mole¬ 
cules  s'1  m“*. 

Some  authors  have  questioned  the  general  validity  of  (I)  on  various  grounds 
[4-7].  However,  in  the  literature  (1)  is  almost  always  successfully  fitted  to  experi¬ 
mental  data.  Indeed,  this  empirical  relation  can  be  given  mechanismic  jutification 
(see  [3]). 

In  dynamic  experiments,  also,  it  is  commonly  agreed  that  (1 )  may  be  adapted  to 
describe  the  reaction  rate: 

i  =  *,/<!- x)  exp  [-£/*  7V)]  (2) 

T(r)  is  controllable  by  the  experimenter.  Some  temperature  programs  offer  the 
mathematical  advantage  that  exp  [-EikT\/t  can  be  integrated  analytically  (see 
later  discussion  on  integral  methods  of  data  analysis).  Examples  are  the  hyperbolic 
program  where  l/T  **  A  —  Bt.  i.e.  t  =  BT'  [8,  9]:  a  parabolic  program  in  which 
AT~  +  BT  -  C  =  /.  i.e.  1  it  =  2 AT  +  B.  with  B  =  AE/K,  by  iteration  [10]:  and 
an  exponential  program  so  that  t  =  exp  { -  BIT)  with  B  i  E  k  [10].  For  the  sake 
of  experimental  convenience,  however,  the  arrangement  is  usually  that  t  *  4> 
(A.  C  and  <t>  represent  constants  in  a  particular  run  of  experiment). 
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However,  it  has  been  taken  by  some  authors  who  object  to  (2).  that 
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where  =  (a)i,0(hern,al.  The  argument  is  then  that  (2)  is  seen  to  be  inade- 

'  *  1 '  r.0 

quate  even  in  the  case  of  d<p  =  0.  since  the  second  term  on  the  R.H.S.  of  (3)  is  non¬ 
zero  but  left  out  in  (2):  see  e  g.  (I  I.  I2|.  It  has  further  been  proposed  [13].  by  a 
derivation  starting  from  (3).  that  (2)  is  correct  only  if  it  includes  the  extra  factor 

[1  +  (l  -  TJT)E/kT] 

in  which  T„  s  T(0).  But  wc  hold  that  (3)  is  unsound.  Given  i.  T  and  </>.  t  is  not 
uniquely  determined  and  therefore  not  a  function  of  these  system  variables:  d*da  > 
>  0.  as  a  cannot  decrease  even  for  negative  d T  and  d$!  Nevertheless,  the  inexact 
differential  dx  can  be  integrated,  if  the  dynamic  process  can  be  treated  as  the  limit¬ 
ing  case  of  a  series  of  time  intervals,  during  which  the  reaction  proceeds  isother- 
mally  according  to  (1)  but  at  the  end  of  each  of  which  T  is  altered,  in  a  time  so  short 
that  during  it  the  sample  is  unchanged.  Along  this  path  P  the  result  is  easily  ob¬ 
tained  [14]: 

i  T 

j  j  je.pc-Mrwr  ,4, 

0  'f  T. 

where  K  =  Kr  exp(-£/AT).  We  must  emphasize  that  (4)  is  not  logically  self- 
evident.  as  is  sometimes  implied  [14]  or  argued  by  mathematical  operations  based 
on  the  presumption  that  a  =  x(T.  t)  [15],  Rather,  it  comes  from  the  assumption 
that  the  reaction  under  study  involves  no  slow  processes,  so  that  x  depends  only 
on  the  present  values  of  a  and  T  (of  is  a  function  of  state),  but  not  on  the  history  of 
the  reacting  system  (c.f.  [16]).  Only  by  this  assumption  (absence  of  memory  effects) 
can  the  dynamic  process  be  treated  as  P.  (Experimentally,  a  temperature  program 
with  temperature  jumps,  which  approximates  P ,  has  been  realized  on  a  thermo¬ 
balance  interactively  controlled  by  computer  [17].)  Equations  (2)  and  (4)  are  of 
course  equivalent.  Their  validity  has  also  been  shown  by  rational’  thermodynamic 
arguments,  in  which  the  functional  relation  x  =  g  (/.  Rt)  is  regarded  as  the  ‘con¬ 
stitutive  equation’  characterizing  the  reaction  system  [18].  In  a  new  direction. 

the  possibility  has  been  suggested  [64]  that  solid-state  reactions  may  be  studied 
by  far  from  equilibrium  thermodinamics. 

On  the  other  hand,  experience  shows  that  apparently  (4)  is  not  always  followed 
exactly.  Consider  a  reaction  being  investigated  by  a  series  of  experiments  conducted 
at  different  heating  rates  4>  but  with  the  same  initial  temperature  T„.  In  (4)  we  see 
that  the  R.H.S.,  for  a  given  upper  temperature  limit  T.  is  directly  proportional  to 
1/^-  Plots  of  the  L.H.S.  vs.  T  should  therefore  all  have  the  same  shape.  It  may  hap¬ 
pen.  however,  that  increasing  departure  from  isomorphism  is  seen  when  experimen¬ 
tal  data  obtained  at  higher  <p  are  so  analyzed.  The  most  probable  explanation  is  that 
the  temperature  change  is  too  fast,  causing  the  temperature  distribution  in  the 
sample  to  become  significantly  non-uniform.  In  fact,  thermal  equilibrium  is  an 
underlying  assumption  when  (4)  is  derived  above:  without  it  i  will  depend  on  the 
thermal  history  of  the  reacting  system. 

Other  factors  may  also  be  at  work.  The  reaction  rate  may  be  sensitive  to  the 
structure  of  the  reactant,  and  a  higher  0  can  enhance  the  defect  density  or  change 
the  activation  energy  of  reaction  at  a  defect  site  [19].  In  branched-cbain  reactions, 
the  speed  of  the  progressive  accumulation  of  active  centres  may  vary  with  the  rate 
of  change  in  temperature  [20].  If  the  reaction  is  a  surface  process,  the  distribution 
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of  reaction  ‘centres'  among  corners,  edges  or  faces  of  the  sample  may  change  with 
<f>  [21].  It  may  also  be  that  the  chemical  system  under  study  has  multiple  reactions 
proceeding  concurrently  in  it,  and  they  have  different  £  [22].  All  these  variations 
in  £  may  be  accompanied  by  changes  of  Kx  in  the  same  direction.  Because  of  this 
coupling,  a  linear  relation  between  £  and  log  AT,  is  sometimes  observed.  Called  the 
'compensation  efTect’,  this  phenomena  does  not  necessarily  mean,  as  was  suggested 
[4],  that  the  Arrhenius  expression  in  the  R.H.S.  of  (4)  is  invalid.  In  all  these  cases, 
by  varying  <f>  the  experimenter  can.  in  fact,  gain  additional  insights  into  the  mecha¬ 
nism  of  the  reaction,  or  distinguish  between  the  competitive  reactions  in  the  react¬ 
ing  systems  (a  situation  usually,  though  not  always,  indicated  by  the  presence  of 
multiple  peaks  in  the  x  curves).  This  is  possible  if  the  method  of  data  analysis 
employed  is  such  that  K„  and  £  are  determined  from  a  single  a  ot  a  curve,  rather 
than  from  data  at  a  number  of  heating  rales.  The  method  wc  suggest  will  be  of  this 
type. 

It  may  also  happen  that  in  (2)  Kx  x  T,  so  that  it  cannot  be  taken  outside  the 
integral  sign  in  (4).  Indeed,  modifications  have  been  suggested  of  some  methods  of 
data  analysis  (those  that  assume  a  reaction  order  for  the  reaction)  to  take  this  extra 
temperature  dependence  into  consideration  [23].  However,  even  when  theoretically 
required,  the  correction  may  for  practical  purposes  be  ignored,  unless  £  is  small  or 
temperatures  used  are  very  high  (d  In  K/dT  «*  [£  +  krykT*).  Likewise,  any  slight 
temperature  dependence  of  £  can  usually  be  neglected.  Furthermore,  irrespective 
of  this  or  the  above  complications  the  form  of/f  1  -  a)  in  (2)  and  (4)  is  not  affected. 
It  should  be  the  same  as  that  in  (I).  on  the  basis  that  the  dynamic  process  can  be 
treated  as  the  limiting  case  of  a  series  of  isothermal  intervals,  as  already  mentioned 
above. 


Data  analysis 

Experimentally,  A  or  a  is  obtained  by  DSC.  DTA,  TG,  DTG,  quantitative  IR 
spectroscopy  or  X-ray  diffraction,  dilatometry,  or  measurements  of  chcmi-lumincs- 
cence.  ultrasonic  attenuation,  dielectric  constant,  viscoelasticity,  thermal  or  electri¬ 
cal  conductivities,  or  optical  reflectivity  when  changes  in  these  characteristics  can 
be  correlated  with  x.  The  oldest  technique  is  thermomanometry,  in  which  the  pres¬ 
sure  of  an  evolved  gas  is  measured  at  constant  volume,  but  its  use  has  so  far  been 
more  popular  in  isothermal  experiments.  Many  methods  of  analyzing  d  or  a  data 
have  been  proposed  to  calculate  the  kinetic  constants  £  and  Km  (for  a  critical 
review  of  the  earlier  work  see  [24])  and  sometimes  also /(I  -a).  Often  they  were 
originally  formulated  with  reference  to  one  particular  instrumentation,  but  they 
may  be  made  generally  applicable  to  all  techniques  after  quantities  measured  on 
DSC,  DTA.  TG  instruments  etc.  are  all  interpreted  in  terms  of  a  and  A.  On  the 
other  hand,  their  validity  does  depend  on  the  particular  reaction  whose  data  are 
being  analyzed.  Their  limitations  in  this  respect  form  the  subject  of  our  discussion 
below.  They  will  be  examined  in  three  groups:  peak-temperature,  integral  and 
derivative  methods,  in  this  order.  Sophisticated  instrumentation  systems  are  com¬ 
ing  into  use.  that  incorporate  computers  to  establish  baselines  or  other  null  settings, 
to  carry  out  automatic  data  acquisition,  and  to  let  the  experimenter  interactively 
analyse  the  data  (e.g.  [25]).  Such  advances  do  not.  however,  remove  the  danger  of 
uncritical  choices  of  the  method  of  data  reduction. 
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Critical  examination  of  current  methods 

Peak-temperature  method 

Kissinger  [26]  considers  reactions  of  the  type /{I  -a)  *  (1  -  a)".  Differentiating 
(2)  with  respect  to  /,  and  setting  the  resulting  expression  to  zero,  he  obtains 

iJE/kUITi)  -  exp(-£/Arrm)/Kl (5) 

in  which  m.  signifies  'p:ak'  quantities,  at  the  point  of  maximum  £  where  a  =  0. 
He  next  a»>umes  that  n(l  —  *mV  1  a  I :  therefore 

4>Tj CL  CKP(- ElkTJ  (6) 

regardless  of  n ,  which  itself  may  be  calculated  from  the  shape  of  the  a  (/)  curve. 
E  ami  k.  on  the  other  hand,  are  obtained  by  performing  a  series  of  experiments  at 
different  <f>.  An  aspect,  which  we  regard  as  an  inefficiency,  of  Kissinger's  method  is 
that  only  one  point  on  the  curve  is  used  although,  in  the  case  where  multiple  peaks 
occur  signifying  that  different /{ I  —  a)  and  E  govern  different  sections  of  the  curve, 
the  method  should  still  be  applicable  to  each  peak. 

There  is,  however,  an  important  limitation.  The  a  priori  condition  that/f  1  -  a)  * 
=  (l-a)"  is  actually  valid  only  in  very  special  circumstances,  namely  when  the 
rate-limiting  step  of  the  reaction  is  the  inward  movement  at  a  constant  speed  of 
the  reactant-product  interface,  where  n  is  0.  >/-.  or  2/3  for  one-,  two-,  or  three- 
dimensional  movement,  respectively,  or  when  the  reaction  is  unimolecular  so  that 
n  —  I.  Even  among  these  special  cases,  the  other  approximation  that  Kissinger 
uses  is  still  conditional,  since  n(l  - 1  s  1  only  for  n  =  1 .  When  n  is  l/2  or  2/3. 
this  expression  varies  with  approximately  as  dam  n(l -n)/(l -am)J  "  ^  0.2 
da„.  where  Aim  is  the  variation  in  <xm  itself.  In  the  Appendix  we  show  that  am 
changes  with  </>  in  the  general  case.  Hence,  when  an  apparent  reaction  order  exists 
and  is  >/>  or  2/3,  Kissinger's  method  can  lead  to  a  systematic  deviation  in  (6)  and  thus 
generate  a  significant  but  hidden  error  in  the  calculated  Eand  K, . 

If  no  apparent  reaction  order  exists,  then  it  definitely  should  not  be  used,  other¬ 
wise  an  approximately  linear  plot  from  (6)  results  in  totally  misleading  values  of 
the  kinetic  constants.  An  example  is  in  the  decomposition  of  benzenediazonium 
chloride:  it  derives  from  DTA  data  a  value  of  E  that  is  40%  lower  than  the  nearly 
identical  values,  obtained  by  applying  other  methods  of  analysis  to  the  data  from 
DTA  as  well  as  other  techniques  [27  J.  Other  examples  are  in  the  study  of  lithium 
aluminium  hydride,  where  the  Kissinger  values  are  half  of  the  isothermal  result 
[28],  in  RDX  where  it  is  again  40%  lower  than  all  the  values  calculated  by  other 
methods  [29],  and  in  urea  nitrate,  where  it  is  30%  lower  [30]. 

Integral  methods 

The  L.H.S.  of  (4)  is  a  function  of  a  only  and  will  be  denoted  by  F(a);  the  R.H.S. 

r 

can  for  practical  purposes  be  equated  with  \  K/<t>  d  T,  since  in  experiments  Tn  will 

b 

be  such  that  reaction  velocity  is  negligible  below  it,  i.e.  T„  <  E/k.  In  view  of  these 
considerations,  many  authors  have  proposed  different  methods  of  analysing  a (T) 
data.  T 

The  temperature  integral  (  exp  (—E/kT)dT  has  no  analytical  solution.  (In  the 

o 

unusual  case  of  a  hyperbolic,  parabolic  or  exponential  temperature  program, 
on  the  other  hand,  exp  (-E/kT)lt  is  integrable.)The  numerical  values  of  the  inte¬ 
gral  have  been  compiled  but,  being  a  function  of  both  E  and  T.  are  not  directly 
useful  unless  an  iterative  solution  of  (4)  by  trial-and-error  is  resorted  to.  Such  an 

approach  has  been  advocated  by  Zsakd  [31]  who  considers  in  particular  the  cases 

* 

of. /II  -at)  =  (I  —  at)*  with  n  —  0,  1/3,  1/2,  2/3,  I  or  2,  when  $(at)«  log  I  ■■ 
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has  simple  analytical  expressions,  and  by  Satava  and  SkvAra,  [32)  who  generalize 
the  method  slightly  by  tabulating  the  values  of  /Ha).  0  <  a  <  1,  for  some  other 
forms  of/)  1  -  a). 

For  more  efficient  approaches,  approximations  to  the  integral  are  necessary. 
Thus,  taking  the  first  two  terms  in  an  asymptotic  (u  s  £/AT  -»  oo)  series  of 

k 

r  <jt 

|  exp  ( —  t)  Coats  and  Redfern  obtain  the  linearized  relation  [33): 

\n(F(z)T-*)  =  A  -  E/k(T-')  (7) 

where  A  =  ln(A  ,  A /£</>)(  I  —  Ik  T/t )  is  ‘sensibly  constant'  if  the  range  of  temperature 
AT  is  small.  They  further  assume  that /(I  -  a)  -  ( I  -  a)*,  and  so  F(a)  can  be  cal¬ 
culated  at  each  (a.T).  Plotting  (7)  for  several  values  of  a  thus  gives  E  and  K 
making  use  of  only  one  set  of  data  corresponding  to  a  single  <)>. 

Several  cautionary  notes  should  again  be  made  here.  The  assumption  for/)  I  -  a) 
has  already  been  discussed.  Similar  to  the  case  of  Kissinger's  method,  results  ob¬ 
tained  may  be  wrong  and  misleading  if  this  functional  form  is  not  independently 
determined  beforehand.  Thus,  in  a  study  on  the  dehydroxylation  of  kaolinite  (34), 
straight  lines  over  different  ranges  of  (T~ *)  are  given  by  (7)  for  a  whole  series  of 
values  of  n.  namely,  0,  0.5.  0.667,  a  and  2.  In  particular,  plots  using  n  *  I  and 
n  -  2  are  almost  equally ‘good’. 

Secondly,  the  accuracy  of  the  asymptotic  approximation  is  rather  low.  By  com¬ 
paring  its  values  with  tabulated  values  of  the  integral  [31.  35.  36),  we  find  its  rela¬ 
tive  errors  to  be  4///  -  20  %  at  u  -  5,  5%  at  it  -  10.  and  1 .5  %  at  u  -  20.  Thus,  for 
example,  if  £  is  I  eV,  then  for  an  accuracy  of  98  %  the  highest  temperature  reached 
in  the  experimental  run  should  not  be  more  than  600  K.  a  very  low  figure  for  most 
materials  though  it  is  higher  for  larger  £.  Additionally,  expanding  A  into  a  power 
scries  shows  that  AA/Asi2k  AT/E, so  that  at  say  2%  inaccuracy  the  range  of  tempera¬ 
ture.  AT,  from  which  (a,  T)  points  are  selected  should  be  less  than  I0OK  (for  £  « 

I  eV).  The  total  possible  deviations  in  the  calculated  £  and  Kr  are.  to  first  approxi¬ 
mation.  the  sum  of  the  All  I  and  A  A/ A.  It  certainly  is  unsatisfactory  if  they  are  largtf 
and  yet  nowhere  mentioned  in  the  calculation. 

Other  approximations  to  the  temperature  integral  have  been  suggested  by  van 
Krevelen  el  at.  (37)  and  by  Horowitz  and  Metzger  [38).  who  made  use  of  certain 
asymptotic  expansions  in  the  vicinity  of  Tm,  the  temperature  at  peak  reaction  rate. 
Both  have  been  shown  [39)  to  be  even  less  accurate  than  the  Coats  and  Redfern 
approach,  and  so  will  be  left  out  in  our  discussion. 

Amongst  the  integral  methods,  the  best  is  probably  the  one  due  to  Ozawa,  which 
requires  data  at  different  4  but,  in  it  /(l  -  a)  remains  completely  general.  The  ap¬ 
proximation  to  the  temperature  integral  is:  - 

r 

J  exp< -ElkTytT  *  /  I  O'  -*.»-*«***n  (8) 

so  that  from  (4) 

log p,  +  0.457(£/A)/r,  -  log*,  +  0.457(£/A)r,  (9) 

where  7",  and  Tt  are  taken  at  an  arbitrary  but  identical  value  of  a  in  the  two  curves 
corresponding  to  heating  rates  di  and  #t*  Plotting  log  #  vs.  1  IT  for  selected  values 
of  a  should  therefore  produce  straight  lines,  the  slopes  of  which  give  £  [40]. 

Three  comments  are  appropriate  here.  By  comparing  (8)  with  tabulated  numer¬ 
ical  values,  we  see  that  it  is  7%  out  at  u  -  10 or  T  -  1I70K,  and  3%  and  less  only 
for  T  <  720K  (if  E  -  1  eV).  These  errors  should  be  examined  before  Ozawa’s 
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method  is  applied.  Secondly,  the  method  has  been  modified  [41]  to  read,  in  place 
of  (9). 

A  In  <t>/ATm 1  -  0.457  Elk  (10) 

in  which  d  I nt)>  s  ln^1  —  In  <j>s,  etc.,  and  m  denotes,  as  before,  peak  quantities. 
This  relation  may  be  compared  with  (6)  but  in  general  it  does  not  hold  since,  as 
shown  in  the  Appendix.  im  varies  with  4>.  Lastly,  like  Kissinger's  method.  E  cannot 
be  determined  from  data  at  a  single  <t>.  and  in  some  cases  this  may  be  adisadvantage. 
as  discussed  before. 

Derivative  methods 

The  derivative  methods  offer  an  advantage  over  those  described  above  in  invok¬ 
ing  no  mathematical  approximations.  Unfortunately,  they  use  &  data  which,  with 
present  instrumentation,  tend  to  be  of  lower  quality  whether  they  are  obtained  by 
numerically  dilftfrcntiating  (he  *  data  or  are  direct  experimental  read-outs. 

The  most  straightforward,  but  as  it  stands  relatively  inefficient,  of  the  derivative 
methods  is  to  write  (2)  as:  — 

In  (a//(  I  -  a))  =  In  K,  -  E/kT  (11) 

and  to  substitute  different  of /(I  -*)  until  a  linear  plot  appears  [42].  Later,  we  shall 
argue,  however,  that  even  this  labor  omnia  vincit  approach  like  all  dynamic  meth¬ 
ods  in  general,  cannot  guarantee  correct  values  of  E  and  K,  (nor  an  unambig¬ 
uous  form  of  /( 1  -  a)  in  this  specific  case),  although  the  labour  it  involves  may  be 
undertaken  by  the  computer. 

The  earliest  derivative  method  is  probably  that  of  Borchardt  and  Daniels,  origi¬ 
nally  formulated  for  homogeneous  reactions  in  the  liquid  phase  [43]  but  later  extend¬ 
ed  to  solid-state  reactions  [44]  for  which  it  is  now  frequently  used.  The  method 
puts /( 1  —  a)  *  ( 1  -  a)"  into  (II),  with  n  given  a  guessed  value,  and  if  a  linear  plot 
results  then  E  and  Ky  are  obtained  from  it.  Based  on  this  method.  Hauser  and 
Field  [45]  have  developed  a  computer  procedure,  in  which  plots  are  generated  for 
a  series  of  values  of  n  incremented  at  discrete  steps,  and  the  ‘best'  one  is  then  se¬ 
lected  to  yield  E.K and  n.  An  attraction  of  this  method  is  that  n  can  be  readily 
selected  by  eye.  Alternatively,  since  in  this  case 

A  In  a  dln(l  -  a)  »  — (£  k[AT~ 1  A  In  (I  -a)]')  +  «  (12) 

a  plot  of  the  L.H.S.  vs.  the  quantity  in  the  square  brackets  at  once  gives  E  from  the 
slope  and  n  as  the  r-intercept[47).  If  constant  din*.  dln(l  -  a),  or  AT~  1  is  selected. 
Eq.  (12)  can  be  further  simplified  [61],  We  have  emphasized  previously  the  falli¬ 
bility  in  presuming  such  a  convenient  form  of/]l  -a):  Ozawa  [47]  has  commented 
on  the  possibility  that  this  procedure,  and  the  integral  method  of  Coats  and  Red- 
fern.  may  give  false  values  of  E  and  K, .  In  addition,  since  (12)  involves  the  ratios 
of  differences,  the  quality  of  data  called  for  is  even  higher  than  that  demanded 
alone  by  the  use  of  a :  experimental  data  so  plotted  more  often  than  not  show  very 
large  scatter.  The  Rogers  and  Morris  method  [48]  plots  din*  i  s.  T~ and  can  be 

seen  to  be  the  special  case  of  n  =  0  in  (12).  An  example  of  the  general  danger  that 
very  linear  plots  may  sometimes  appear  even  if  the  applied  method  is  not  valid  is 
given  by  Patel  and  Chaudhri.  The  Rogers  ar  d  Morris  method  was  used  to  analyze 
DSC  data  on  lead  azide,  and  a  straight  line  results  although  the  calculated  E  turns 
out  to  be  180°o  larger  than  the  Ozawa  value  [49].  Conversely,  the  coincidence  of 
values  calculated  by  various  methods  need  not  prove  that  these  methods  8re  all 
applicable  to  the  case  in  hand.  A  counter-example  is  provided  by  a  DSC  study  on 
RDX  [29],  where  the  Rogers  and  Morris  value  agrees  well  with  other  valuesbut 
the  complex  decomposition  is  beyond  doubt  far  from  the  n  »  0  type. 

On  the  other  hand.  DAvid  and  Zelenyanszki  [50]  plot  In  j  ^  (I  —  at) j/(  1  -a) 

against  (7"  ’);  this  amounts  to  assuming  a  ‘reaction  order'  n  =  I.  It  serves  as  yet 
another  example  of  the  futility  of  linear  plots,  for  their  method  gives  such  plots  for 
the  decomposition  of  'a  wide  range  of  materials'  including  calcium  oxalate  and 
polyethylene  which,  most  likely,  are  not  of  first  or  any  other  'order'. 
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Some  of  the  integral  and  derivative  methods  described  in  the  foregoing  have  been 
compared  by  testing  their  accuracies  on  synthesised  DTA  data  (exact  as  well  as 
with  artificial  random  error)  for  one  E  value  and  temperature  range,  the  reaction 
considered  being  of  the  type  with  a  reaction  order  [51].  Among  the  methods  not 
included  there  is  thatdueto  Friedmann  [52].  It  probably  is  the  most  general  among 
the  derivative  methods.  Like  Ozawa's  procedure,  it  makes  no  assumption  about 
/( I  -  a),  although  it  requires  A  data  which,  furthermore,  have  to  beat  a  number  of  </>. 
Once  again,  from  (2)  with  da  d T  =  xr: 

In  (aTd>)  =  In  (Kxf(  I  —  a))  -  EjkT.  (13) 

Since  K  ,  )\  I  —  a)  is  identical  for  the  same  value  of  a.  taking  a T  and  the  correspond¬ 
ing  T from  several  tj>  one  can  determine  E. 

It  is  our  contention  that  even  Friedmann's  method  has  one  basic  limitation  which, 
more  significantly,  is  shared  by  all  dynamic  methods  described  above.  The  point 
in  question  is  that  all  of  them  have  to  presume  the  constancy  of  f(  \  -  a)  as  the  tem¬ 
perature  is  changed.  However,  since  mechanisms  of  solid-state  reactions  are  gener¬ 
ally  complicated,  there  is  no  general  justification  for  this  presumption,  though  it 
may  be  true  for  particular  reactions  within  specific  temperature  ranges.  An  illustra¬ 
tion  is  the  case  where  parallel  reaction  paths  exist,  each  with  values  of  Kx  and  E 
such  that  a  quantitative  change  of  T  will  lead  to  a  qualitative  change  in  the  dominat¬ 
ing  path.  Another  case  is  where  the  identification  of  the  rate-limiting  step  depends 
on  T.  Methods  have  been  proposed  which,  by  the  use  of  computers,  try  different 
forms  of  /(I  -  a)  in  analyzing  the  dynamic  data  [53-55].  However,  the  search  is 
limited  to  functional  forms  which  are  already  known. 

More  importantly,  from  our  own  experience  with  azides  we  have  strong  doubts 
as  to  the  exactness  in  determining/(  1  —  a)  or  even  its  constancy  from  dynamic  data. 
Likewise,  in  a  study  on  the  dehydration  of  manganese  formate  [54]  for  instance, 
no  unique  form  of  /( 1  -  a)  and  correspondingly  no  unique  values  off  are  identified 
even  over  appropriately  restricted  ranges  of  a,  the  criterion  used  being  minimum 
standard  deviation  in  the  Arrhenius  plot.  Further  examples  are  the  thermal  dehy- 
droxylations  of  kaolinite  [34]  and  of  magnesium  hydroxide  [56],  We  suggest  that, 
in  dynamic  experiments  since  data  are  collected  under  variable  temperature  condi¬ 
tions,  the  change  due  to  /(I  — «)  is  inherently  masked  by  that  due  to  AT(T).This  pitfall 
is  illustrated  in  the  Figures.  Figure  la  shows  the  graphs  of  x{T)  and  its  derivative 
which  are  generated  artifically  according  to  the  theoretical  equation  »T  *  Kx 


10J  IT,  K" 


Fig.  ta.  Artificial  data  t  =  Kt,  and  corresponding  i'  e  dr/dT data,  plotted  against  T  which 

rises  linearly  with  f 

Fig.  lb.  Arrhenius  plots  of.  (1)  #r/3(l  -  i)'T  (2)  *r/2<  1  -  i)1'*.  (3)  >r.  and  [(4)  jr/3j*/* 

for  the  data  shown  in  Fig.  la 
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exp  <  -EjkT)  throughout  from  a  =  0  to  a  =  1  •  The  values  chosen  for  Kx  and  E  are 
108  and  leVrespectively.  Let  us  now  examine  how  the  data  generated  according  to 
this  relation,  which  is  of  the  type  a  =  Kt ,  will  be  fitted  by  different  kinetic  equations, 
one  of  them  being  the  correct  one.  In  Fig.  lb  we  plot  against  lOr'/T  the  natural 
logarithms  of  the  following  expressions:  — 

(1)  ar/3(l  —  a)2'3,  i.e.  assuming  1  -  ( I  - a)J,;‘  =  Kt:  reaction  controlled  by 
three-dimensional  contraction  of  phase  boundary : 

(2)  ar/2(l -a)'  *,  i.e.  1  -  (1-a)'  *  =  Kt  type: 

(3)  ar:  the  original  assumption :  and 

(4)  y'/3ot~':\  i.e.  a  =  ( Kt f :  reaction  controlled  by  e.g.  three-dimensional 
growth  of  existing  nuclei. 

It  is  seen  that  the  incorrect  f\ I  —  a)  in  ( I )  and  (2)  still  give  virtually  linear  plots,  with 
slightly  different  slopes:  interestingly,  curve  (4)  is  so  misleading  as  to  show  two 
“linear”  segments  with  a  seemingly  significant  transition  in  between.  Experimen¬ 
tally,  Guarini  etal.  [57]  have  noted  that  it  is  impossible  to  ascertain  from  their  DSC 
data  whether  the  monomerization  of  9-Me-10-AcAD  has  an  apparent  reaction 
order  of  1.  0.67,  or  0.5,  in  all  of  which  cases  £  has  about  the  same  derived  value. 

The  suggested  approach 

In  the  foregoing  sections,  we  have  discussed  the  limitations  regarding  the  applica¬ 
bility  of  various  methods  that  have  been  used  to  analyse  dynamic  data.  In  many 
published  works  we  find  that  often  a  number  of  apparently  different  methods  are 
used  to  analyse  the  same  set  of  data.  However,  we  think  that  in  many  cases  this 
procedure  is  of  no  real  significance,  when  some  of  the  methods  used  are  mathemat¬ 
ically  equivalent  and  therefore  lead  to  the  same  results,  or  when  some  are  invalid 
in  the  given  situation  and  thus  lead  to  doubtful  values.  The  limitations  of  the  meth¬ 
ods  express  themselves  both  as  discrepancies  in  the  calculated  values  of  the  kinetic 
constants,  and  sometimes  as  fortuitous  agreements  when  some  of  the  methods  are 
certainly  inapplicable.  (An  extreme  example  of  the  second  situation  is  that,  for 
RDX,  the  Kissinger  value  [29]  of  £  is  near  to  that  obtained  [58]  by  plotting 
7  rs  T~\  a  procedure  which  has  absolutely  no  theoretical  justification.)  Accord¬ 
ingly.  we  suggest  that  the  interpretation  of  dynamic  data  should  as  far  as  possible 
be  based  on  results  from  isothermal  experiments.  A  similar  approach  has  been 
used  for  studying  the  dehydroxylation  of  kaolinite  by  Achar.  Brindley  and  Sharp 
[63]. 

One  can  unambiguously  determine/!  1  -  a)  over  the  whole  range  of  a  and  over 
the  relevant  temperature  range,  from  the  independent  analysis  of  individual  iso¬ 
thermal  curves.  A  systematic  method  of  efficiently  implementing  this  identification 
has  been  proposed  by  us  [3].  It  may  also  be  noted  that  thermoanalytical  equipments 
are  equally  applicable  in  isothermal  experiments  (see  e.g.  [62])  though  they  are 
more  often  used  in  the  dynamic  mode.  The  identified  form(s)  of/[l  -  a)  can  then 
be  substituted  into  either  (2)  or  (4).  In  this  way.  form  the  dynamic  a  or  &  data  one 
can  then  determine  accurately  the  non-average  and  single-sample  values  of  £  and 
Kr :  advantages  which  have  been  mentioned  in  the  introduction  to  this  paper. 
Moreover,  the  values  will  correspond  individually  to  different  heating  rates. 

We  applied  this  approach  to  the  spinel  formation  ZnO  +  CrtO,  -» ZnCr,0«. 
A  DTA  curve  (experimental  atmosphere:  N,  at  300  mm  mercury)  was  published 
in  Ishii  et  at.  [59],  who  have  also  monitored  a(0  by  chemical  analysis  when  the 
reaction  proceeded  isothermally  in  nitrogen  flowing  at  50  ml/min,  and  showed 
that  the  isothermal  data  fit  [1  -  (1  -x)l,3]J  -  Kt.  We  have  measured  K  from  the 
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Isothermal  data 

r. 

deg.  C 

A 

rctetive  units 

In  A 

800 

1 

0 

900 

4 

1.4 

1000 

23 

3.1 

Tabic  2 

Dynamic  data 

/. 

,  > 

In  I  X. t  II 

deg.  V 

relative  units 

t  II  a) 

700 

0.1  2 

—  2.6 

820 

0.15  4 

-1.4 

900 

0.2  7 

-0.54 

1000 

0.35  12 

+0.76 

experimental  data  points  at  t  =  20  min  in  the  published  isothermal  plots.  From 
these  values  of  K,  given  in  Table  1.  we  calculate  a  value  of  1.5  eV  for  £. 

In  Table  2  the  values  i  were  measured  from  the  published  DTA  curve  whose 
heating  rate  was  unspecified,  and  the  a  values  were  read  off  from  the  (a,  T)  graph 
which  Ishii  el  at.  have  drawn  presumably  by  integration.  Now,  from  their  analysis 
of  the  isothermal  data  the  governing  kinetic  equation  is,  in  differential  form,  &  — 
=  Af/[(l  -a)"*'*  -  (1  -3t)',,s).  at  least  within  the  ranges  800- 1000°  and  a  from 
0  to  se  0.6  corresponding  to  Kt  *  0  to  &  0.07.  The  Arrhenius  plot  of  a  (1  -  a)"  1,5 
Rl 1 1  for  the  four  data  points  shown  in  Table  2,  is  indeed  a  good 
straight  line.  From  the  plot  we  obtain  £  *  1.3  eV.  In  view  of  the  probable  experi¬ 
mental  errors  and  inaccuracies  in  obtaining  data  from  the  published  graphs,  we  con¬ 
sider  satisfactory  the  reasonable  agreement  between  this  value  and  the  one  calculat¬ 
ed  from  the  isothermal  data. 


Conclusion 

Most  of  the  commonly  used  methods  of  analyzing  dynamic  data  have  been  shown 
to  be  applicable  only  under  particular  conditions.  It  has  been  pointed  out  that  to 
use  these  methods  without  considerations  of  the  range  of  their  validity  can  give 
misleading  values  of  the  kinetic  parameters.  An  approach  has  been  advocated  in 
which  use  is  made  of  both  the  dynamic  and  isothermal  data:  the  functional  form 
/(l  -2)  is  determined  from  the  isothermal  experiments.  This  form  in  conjunction 
with  the  dynamic  data,  gives  the  values  of  the  kinetic  constants. 
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Appendix 

We  wish  to  predict  how  the  value  of  x  at  peak  reaction  rate  varies  with  the  heat¬ 
ing  rate  0.  At  a  =  am:  from  (4) 


da 

J  AT^) 


0  T. 

and  from  the  fact  that  a  =  0 


J  dTexp(-ElkT)=  F(*m.  0.  TJ  =  0 


rn  ,  .  _  E  $  cxp(  E/kTm) 

'  ^  ~  ^a'"'  <t>'  t)  —  0  • 

Solving  the  simultaneous  equations  dF  =  0  and  dC  =  0.  we  find 


(14> 

(15) 


da„,  /  t'F  i)G  iF 

t'G  |  /  ,'F  iXi  i'F  i(,  ) 

d 0  =  |<>  <‘7m  t  T„, 

•'4>  '  1  <'Tm  ia„,  i)T„  j  ’ 

(16) 

Defining  dimensionless  quantities 

V  =  ElkTn.  n  s  <j>iK  ,Tm.  and 

^  \ 

da //( 1  —  a),  we  have  the  following : 

0 

cF  I 

cG 

/d-aj  ' 

- - /<  1-0. 

?F 


cr  a.  _ 
r,v  =  -t‘-  dTcxp(-U)  = 
'  0  0 


cF 


To 

I 


l(y,„)  cG  nl 

“’7#  =  1>  CXP(l/)- 


cg  ftua  +jj) 

~  Tm 


<  Tn  ncxp(U)Tm '  erm 

Hence  —* "  -  -  ^  J*2  +  WJ  -  1  _ 

dd>  <t>  h0( 2  +  (/)//(  1  -  a,„)  -  /*<  I  -  ,„)jh 


exp(U).  (17> 
(18) 


where  h  =  /i  exp  ((/).  Incidentally,  dTJd<t>  can  tc  derived  in  a  similar  way. 

The  only  case  we  find  reported  in  the  literature,  in  which  am  is  apparently  inde¬ 
pendent  of  0,  is  the  primary  recrystallization  cf  pre-cctr.prcsscd  copper  [60],  where 
xm  =  0.5.  In  all  other  cases,  experiments  give  changing  am.  We  have  made  a  rough 
check  on  (18)  by  taking  the  case  of  the  decomposition  of  the  explosive  RDX  [29). 
for  which  the  Rogers  and  Morris  method  eives  E  =  2.10  eV  and  K,  *  l01Mi~ l. 
The  reaction  is  complex,  but  these  representative  values  are  chosen  because  they 
correspond  to  an  assumed  kinetic  equation  in  which  /)  1  -  a)  =  1 .  We  thus  have 
very  simply  /(am)  =  %„  and/'(!  -am)  *=  0.  For  0  =  0.167  K  s "  K  am  is  given  as  0.62 
and  T„  as  512K:  our  calculation  shows  dam/d0  s  -0.1  sK'1,  a  value  which 
compares  well  with  the  experimental  indication  that  da mIA(f>  =  (0.60  -  0.62)/ 
/ (0.333  -  0.1 67) sK'1. 
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